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Abstract 

A new version of the Teukolksy Master Equation, describing any massless field of 
different spin s = 1/2,1,3/2,2 in the Kerr black hole, is presented here in the form 
of a wave equation containing additional curvature terms. These results suggest a re- 
lation between curvature perturbation theory in general relativity and the exact wave 
equations satisfied by the Weyl and the Maxwell tensors, known in the literature as the 
de Rham-Lichnerowicz Laplacian equations. We discuss these Laplacians both in the 
Newman-Penrose formalism and in the Geroch-Held-Penrose variant for an arbitrary 
vacuum spacetime. Perturbative expansion of these wave equations results in a recur- 
sive scheme valid for higher orders. This approach, apart from the obvious implications 
for the gravitational and electromagnetic wave propagation on a curved spacetime, ex- 
plains and extends the results in the literature for perturbative analysis by clarifying 
their true origins in the exact theory. 
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§1. Introduction 

In recent years there has been considerable interest in perturbations of the gravitational 
field and in other spin s massless test fields (s = 1/2, 1,3/2,2) imposed on a gravitational 
background (especially on a black hole background, which is of particular interest for study- 
ing features of gravitational collapse and consequent gravitational wave emission). The 
approaches to perturbations existing in the literature can be divided essentially into two 
classes: 1) metric perturbations and 2) curvature perturbations. 

Pioneering work in the first class was done by Regge, Wheeler and Zerilli, who produced 
important results for spherically symmetric black holes i^'H'l) by decomposing the metric 
perturbations into multipoles which satisfy a set of (linearized) equations@''&'i ) 'B'i' ) '&. How- 
ever, a) perturbing the metric components written with respect to some coordinate system 
is not a "gauge invariant" procedure as pointed out by Stewart and Walker 0) and b) apart 
from spherically symmetric background solutions, this approach does not seem to be very 
useful. 

Working in the second class, i.e., perturbing directly the components of the curvature 
tensor, one can use frames adapted to the null principal directions of the Weyl tensor of the 
background spacetime: this allows special simplifications (even partial decoupling) of the 
equations themselves. However, this approach has problems too: perturbing the components 
of the curvature tensor with respect to some special adapted frame is not a "tetrad invariant" 
procedure and one must add this difficulty to that of the gauge invariance which is still 
present. 

The "state of the art" regarding curvature perturbations is still represented by the work 
of Teukolsky0> (in the context of the Newman-Penrose formalism 0', NP in the following), 
which partially received its mathematical foundation by Stewart and Walker 0) and was 
subjected to some extensions by other authors 0)'O)'il)>lil)>[i3)>i3>>il)>iI>. Very recently, a 
gauge-invariant second-order perturbation study for a type D vacuum geometry was initiated 
by Lousto and Campanelliil', also stimulated by relevant results concerning the second order 
perturbations of a Schwarzschild black 

We introduce a new form of the Teukolksy Master Equation, closer to a "generalized 
Klein Gordon" equation. This allows us to investigate exact wave equations on a general 
vacuum spacetime: the de Rham-Lichnerowicz equation satisfied by the Riemann tensor 
as well as the ordinary de Rham equation satisfied by the Maxwell tensor in a test field 
approximation. All the results are formulated for tensor components with respect to any 
frame and then specialized to the Newman-Penrose formalism and to its further extension, 
the GHP formalism, introduced by Geroch, Held and Penrose©. Starting from the exact 
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equations for the Riemann and Maxwell tensor components, it is easy to give a unified and 
complete picture of the existing theory of perturbations to any order, framing all the known 
results (which are scattered in many different formalisms) in a more general form. In type 
D spacetimes all the known results, including those of Teukolsky0>, Stewart and Walker0> 
and Lousto and Campanelliii* are easily recovered and discussed. Finally, some remarks are 
given about how this approach can be generalized to half-integer spin test fields as well as 
to the case of nonvacuum spacetimes. 



§2. The Newman-Penrose formalism: a brief review 

In this article we introduce the Newman-Penrose form of the generalized de Rham wave 
operator for the Weyl, Riemann and Maxwell tensors. To accomplish this, it is first nec- 
essary to specify all the formalism conventions: we will follow exactly the notations of 
Chandrasekhar 0). For the sake of completeness we recall here some details. 

A Newman-Penrose frame is defined by four complex null vector fields ei = 1, e 2 = n, 
e 3 = m, e 4 = m*, satisfying the relations 



1 • m = 1 ■ m* = n ■ m = n • m* 
l l = n n = m m = m* • m* 
1 ■ n = 1, m • m* = —1, 

so that the components of the metric tensor in this frame are 

fO 1 \ 
1 


Vo o 



0. 
0. 



(2-1) 



Vab = V 



ab 








-1 





-1 

J 



(2-2) 



With an abuse of notation (not distinguishing vectors and 1-forms) the dual frame is 



n = e 2 , 

— m* = — e 4 , 



e z = 1 = ei 

4 

e = — m = 



-e 3 • 



(2-3) 



The basis vectors, thought of as directional derivative operators when acting on scalars, are 
denoted by 



ei = e 2 = D, e 2 = e 1 = A , 
e 3 = — e 4 = 5, e 4 = — e 3 = 5*, 



(2-4) 
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and the 12 complex spin coefficents are introduced through the following linear combinations 
of the 24 Ricci rotation (connection) coefficents 7 ca b = e c • V e( ,e a = —^ ac b 



K 


= 7311, 


P 


= 7314, 


e = 


2l72ii + 734i), 


a 


= 7313, 


P 


= 7243, 


7 = 




A 


= 7244, 


T 


= 7312, 


a = 


7^(7214 + 7344), 


V 


= 7242, 


7Y 


= 7241, 


P = 


2(7213 + 7343) 



As a general rule, complex conjugation of null frame components is equivalent to interchang- 
ing the indices 3 and 4. Using this prescription the inverse relations of (2-5) are 



7121 


= -(c + 0, 


7122 


= -(7 + 7*), 


7123 


= -(P + a*), 


7124 


= -(a + F), 


7131 


= -K, 


7132 


= -T, 


7133 


= -v, 


7134 


= ~P, 


7141 


= -« , 


7142 


= -r*, 


7143 


= ~P\ 


7144 


= -o- , 


7231 


= 7T , 


7232 




7233 


= X\ 


7234 


= P , 


7241 


= 7T, 


7242 




7243 


= P, 


7244 


= A, 


7341 


= (^-0, 


7342 


* 

= 7-7 , 


7343 


= 13 -a*, 


7344 


= a-(3*. 



(2-6) 



The commutation rules (Lie brackets) [e a , e b ] = C c a be c reduce to 

[A D] =(7 + 7 *)D + (e + e*)A - (r* + tt)5 - (r + 7r*)<T, 

[5, £>] =(a* +0- ti*)D + kZ\ - (p* + e - e*)5 - <7<T, 

[5, Z\] = - + (r - a* - /3)A + (// - 7 + 7*)^ + A*5*, ^ ?) 

[5*, 5] = (,/ _ p)D + (p* - p)Z\ + (a - /T)<5 + (/3 - a*)<5*. 

The 10 independent components C abcd of the Weyl tensor are represented by 5 complex 
scalars in the Newman-Penrose formalism 

% = - C 13 i 3 = - C abcd l a m b l c m d , 
*i = - C1213 = - C abcd l a n b l c m d , 

&2 = ~ C1342 = - C abcd l a m b m* c n d , (2-8) 
% = - Ci 242 = - C abcd l a n b m* c n d , 
&4 = ~ C 242 4 = - C abcd n a m* b n c m* d , 
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with the additional properties 



C1334 


— C1231 — 01, 


Cl241 


— C1443 — ?f r *, 


C1212 


=^3434 = + 


Cl234 


= (^2-^*), 


C2443 


= — C*1242 — 03, 


Cl232 


= ^2343 — — 



(2-9) 



and C1314 = C2324 = C1332 = C1442 = 0. Analogously, the 10 independent components R a b of 
the Ricci tensor are represented by the following (four real and three independent complex) 
scalars, packaged in their tracefree {&) and pure trace (A) parts 



^00 = 


- \Rii, 


022 


- ~ lR22, 


^02 = 


1 

- 2^33, 


020 


= - 2^44, 


011 = 


— ^(-^IS + R34), 


001 


= " ^13, 


0io = 


- \Ru, 


012 


= " 2^3, 


021 = 


- -R24, 


A 


24 



(2-10) 



— (^12 - R34) 



The components R a bcd of the Riemann tensor are then related to the Weyl and Ricci tensors 
by: 



-R1212 


= _ ^ 2 _ - 2<P n - 10/1, 


-Rl324 


=0 2 + 2A, 


-R1234 


=02 - 2 *, 


-^3434 


= - & 2 - 2 * + 20ii - 10/1, 


-R1313 


= -0o, 


-^2323 


^4 , 


-R1314 


= - 0oo, 


-^2324 


= - 022, 


-R3132 


=002, 


-Rl213 


= -!^i-3>oi, 


-R1334 


= -#*-$>„!, 


-Rl223 


=0 3 * + 0i2, 


-R2334 


= - 3 * - 012 , 







(2-11) 



with the additional complex conjugate relations obtained by interchanging the indices 3 and 
4. 

The six independent components F ab of the (antisymmetric) Maxwell tensor are repre- 
sented by the three complex scalars 

00 = 013 , 01 = ^(012 + 043) , 02 = 042 , (242) 
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with the inverse relations 



ab 



V 



o 

'A-<k 

-00 

-<t>l 







*5 \ 

02 





-<t>*l + 01 



J i - 9l 




(2-13) 



Finally the Ricci and the Bianchi identities are given explicitly in the book of Chan- 
drasekhar© (section 1.8). Here they will be referred to according to the components of the 
Riemann tensor which give rise to the corresponding equation: [R a bcd] stands for the Ricci 
identity associated with R a bcd and R a b[cd\e} = for the Bianchi identities. 

Exchange of 1 <-> n, m «-> m* in the various Newman-Penrose quantities implies 



D = Pd^ n"^ 



m— m' „ ^ 



< > m ^On 



A, 



5\ 



(2-14) 



l^n m^m* 
K = 7311 < ► 7322 < > 7422 = ~7242 

1^+n m<-^m* 
T = 7312 < > 7321 < > 7421 = ~ 7241 

l^n m^m* 

0" = 7313 < *• 7323 < ► 7424 = "7244 

1 n m^->m* 
P = 7314 < ► 7324 < ► 7423 = "7243 



-7T, 



-A, 



(2-15) 



e = ^(7211 + 734i) ^> ^(7122 + 7342) 



: (7122 + 7432) = -7, 



2 V ' ' y 2 
a = 0(7214 + 7344) ^^-(7124 + 7344) -(7123 + 7433) = -@, 



(2-16) 



#0 = 


— C1313 


<— 


— > 


#1 = 


— C1213 


1- 

<— 


— > 


#2 = 


- C1342 


1- 

<— 


— > 



— ^2323 < ► — L-2424 

— L-2123 < ► — L-2124 



— (^2341 < > 



-a 



2431 



#4, 

— C*1242 = $3j 
— C*1342 = $2 



(2-17) 



TP l—^ P m^m p 

-F13 < ► r-).?. < ► r; 



23 



21 



>2- 
1* 1 



>1 = -(^12 + ^43) " + ^43) ^21 + ^34 ) 



(2-18) 



These operations are useful since they may be used to extend a given component calcu- 
lation to other components. 
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3. NP gravitational perturbations 



The study of perturbations of the various fields in the NP formalism is achieved by 
splitting all the relevant quantities in the form / = l A + l B + #4 = \P A + tyf + a — 
a A + a B + D = D A + D B + etc., where the A terms are the background and the f?'s are 
small perturbations etc. The full set of perturbative equations is obtained inserting first of 
all these splitted quanties in the basic equations of the theory (Ricci and Bianchi identities, 
Maxwell, Dirac, Rarita-Schwinger equations etc..) and keeping only first order terms. After 
some "ad hoc" algebraic manipulations, which generate second order differential relations, 
one usually obtains second order coupled linear PDE for the curvature quantities. However 
for some very special spacetimes, as for Petrov type D ones (i.e. Kerr metric), some PDE's 
are decoulpled. In the following we'll present an example of this procedure, applied to a 
generic vacuum type D metric, for which it follows that: 

y A = & A = & A = V A = (3-1) 

K A = a A = v A = \ A = . (3-2) 

We start from the exact NP relations: 0: 

(S* - Aa + ix)% -(D-Ap- 2e)lPi - 3^ 2 = (5 + tt* - 2a* - 2(3)<P 00 

-(D - 2e - 2p*)<P 01 + 2a<P 10 - 2k<P u - k*$ 02 (3-3) 

(A - 4 7 + p)% -(5 -At- 2/3)^i - 3(T^ 2 = (5 + 2vr* - 2/3)# i 

— (D - 2e + 2e* - p*)<P 02 - \*<P 00 + 2a<P n - 2k<P 12 (3-4) 



(D - p - p* - 3e + e*)a - (5 - r + n* - a* - 3/?)k - *b = . (3-5) 

In the tetrad language, (3-3) represents i2i3[i3|4] = 0, (3-4) represents i?i3[i3|2] = (they 
are two components of the Bianchi identities) and (3-5) is the Ricci identity component 
-R1313 = 0. In this notation we have: 



OO = --R^ FF = AnT^FF = A<kT u (3-6) 



and in the same sense one has to look at the other <Pij. Obviously, in (|3-6|) tt is not a spin 



coefficient but it is the usual mathematical constant coming from the right hand side of 
the Einstein equations. Because we are in vacuo, after the perturbative expansion, in these 
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relations we have that \P A , \P A , a A , n A , A and all the <& A m are zero and consequently we 
get: 

(5* - Aa + n) A V B -(D-Ap- 2e) A V B - 3k b V a = 
An[(6 + vr* - 2a* - 2(3) A T B - (D - 2e - 2p*) A T* ] 



(3-7) 



(A - 4 7 + - (<J - 4r - 2/5) A ^f - 3a B ^ = 

4tt[(5 + 2tt* - 2/3) - (D - 2e + 2e* - p*) A T r B ' 



(3-8) 



mmJ 



(D - p - p* - 3e + e*) A cr B - (5 - r + n* - a* - 3(3) a k b - V B = . (3-9) 

In the following, we'll obmit the A superscript for the background quantities to simplify the 
notation. The type D background metric satisfy the relations: 

D^ 2 = 3p^ 2 (3-10) 

6& 2 = 3t& 2 , (3-11) 
so multiplying (3-9) by \P 2 and taking into account that 

V 2 Da B = D(V 2 a B ) - a B D& 2 = D(V 2 a B ) - a B 3p& 2 (3-12) 
& 2 5 K B = 5(V 2 k b ) - K B m 2 = 5(V 2 k b ) - k b 3t& 2 (3-13) 
we obtain the relation: 

(D — 3e + e* — Ap — p*)V 2 a B - (5 + tt* - a* - 3(3 - At)V 2 k b - V B V 2 = 0. (3-14) 

We have now to eliminate \P B from the eqns (3-7) and (3-8). To accomplish this task we use 
the commuting relation found by Teukolsky and valid for any type D metric: 

[D-(p+ l)e + e* + qp - p*] (5 - p(3 + qr) 

-[6 - (p + l)(3 - a* + tt* + qr](D - pe + qp) = (3-15) 

where p and q are two generic constants. We apply now (D — 3e + e* — 4p — p*) on (3-8) and 
(5 + it* — a* — 3(3 — At) on (3-7) and we subtract the two resulting equations. We can now 
use (3-15) with the choice p = 2 and q = —4 to eliminate the \P± term. The combination o B 
and n B remained is exactly the same present in eq. (3-14), and consequently these quantities 
can be eliminate in favor of the $2$cf term. The final equation is: 

[(D - 3e + e* - Ap - p*)(A + p - 4 7 ) - (5 + vr* - a* - 3(3 - At) 

x (5* + tt - Aa) - W 2 ]$o = 4ttT (3-16) 
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where we have posed: 

T = (6 + it* - a* - 3/3 - 4r) [(D - 2e - 2p*)T l B m - (5 + it* - 2a* - 2(3)T t f] 
+ ( D - 36 + e* - Ap - p*)[(5 + 2n* - 2/3)l£ - (D - 2e + 2e* - p*)T* m ] . 

(3-17) 

In the same way one can build the other relations to perform the perturbative analisys 
of a type D vacuum spacetime. We point out that the same technique can be applied to 
the NP Maxwell equations to decouple them as well to the neutrino and Rarita-Schwinger 
onesS-*'© It is important to stress that this perturbative theory has a purely "algebraic and 
manipulative" nature, without taking into account the geometry underlying the problem. 
In the following we'll give to this subject a new geometrical meaning, framing it into the 
context of the field theory. But to accomplish this task we have to adapt the previous studies 
to the Kerr solution, which describes a rotating black hole of specific angular momentum a 
and mass M, with the condition of flat spacetime at infinity. 



§4. The Teukolsky Master Equation 

The Kerr solution in Boyer-Lindquist coordinates is given by§>: 

2Mr 



E 



, 2 AaMr sin 2 9 1 1 , E 2 
at H — dtdcp — —dr 



Ed9 2 - 2Mr 



r 2 + a 2 + 



a 2 sin 2 9 



sin 2 9d<p 2 



where as usual: 



(4-1) 



A = r 2 - 2Mr + a' 



E = r + a cos 2 9 



(4-2) 



For any Petrov type D metric, and in particular for the Kerr solution, simplifications occur 
working in a tetrad adapted to the two repeated principal null directions of the corresponding 
Weyl tensor. In the NP formalism the following quantities 0* completely describe the Kerr 
solution (in this section we use the label A over some quantities in view of a perturbative 
analysis, as will be clear in the following). The Kinnersley tetrad®: 



A 
1 



[r 2 + a 2 ,A,0,a] 



in 



2E 



[r 2 + a 2 , -A,0,a] 



(4-3) 



V2(r + 



ia cos i 



[ia sin 9, 0, 1, 



sin^ 
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with the 4th tetrad vector being (m*^) A , the complex conjugate of [m^) A . The Weyl tensor 
is represented by: 

= $ A = $r A = ^ = 

(44) 

$ A = M{p A f 

while the Ricci tensor and the curvature scalar are identically zero because of the vacuum 
spacetime. The spin coefficents are given by: 



K A = a A = X A = v A = t A = 



A n *A , 



A —I A -iap A p* A sin 9 

^ (r — ia cos 9) ' ^[2 

nA -p* A cot6 a ta(p A ) 2 sm6 

P = n = V2 ' (4 ' 5) 

, (p A ) 2 p* A A A A p A p* A (r - M) 

V = g ' 7 =/i + 2 ' 



. A „A q*A 



a =ii — p 



It is a relevant result that all the massless (non trivial) perturbations of different spin of 
the Kerr metric are included into the Teukolsky Master Equation (TME) 0) which has the 
following form (jr is the mathematical constant): 



(r 2 + a 2 ) 2 



-A 



-2s 



A 

Or 



a 2 sin 2 9 



d 2 ib AMar d 2 ib 
— - H — + 

dt 2 



a 
A 



d 2 tfj 



A 



s+l 



dr J 



- 



Mir 2 - a 2 ) 



l_d_ 
sin 9 89 

i a cos 9 



sin 9 



2s 



sin 2 9 
a{r - M) 
A 



+ 



A dtd<p 
~d9J 

?± + (s 2 cot 2 9 -s)j> = 4n£T. 
at v J 



9 

i cos 9 
sin 2 9 



dtp 



(4-6) 



The following table shows the various quantities involved in the Teukolsky Master Equa- 
tion (T a and J a are respectively gravitational and electromagnetic sources): 



TME quantities 
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S 


T 




<p 





see ref.l^ 




Xo 


+ 1/2 


see ref.l^ 






- 1/2 


see ref.0 




0o 


+ 1 


Jo 


p 




- 1 








+ 3/2 


see ref.N* 


p" 




-3/2 


see ref.ii 




% 


+ 2 


2T 


p 


"%4 


- 2 


2p~ 4 T 4 



(4-7) 



Table 1: The physical field component tp, spin weight s and source terms T for eq. (4-6) 

It is possible to prove that the other NP quantities not described by this equation are 
related to trivial perturbations in the Kerr background©©, and consequently the study 
of the Teukolsky Master Equation alone allows practically the recovery of all the physical 
quantities. The reader can easily verify that eq. (3-16) coincides with (4-6) for s = 2. It is 
remarkable the separability of this equation in the form: 



ip{t,r,9,cf>) 



-*"*e im *R(r)y(0), 



(4- 



by using the so called "spin weighted spheroidal harmonics" which, as it is well known from 
molecular physics^, satisfy the following spectral problem: 



[H + H 1 )Y{6) = -EY{9) 



(4-9) 



with 



Ho 



1 


d 


- ( 


sin£ 


ue 






Hi 


2 2 2 a 

= a oj cos U - 



vn? + s 2 + 2ms cos 9 
sin 2 # 



(4-10) 



(4-11) 



and E = E a j, >m and Y = Y s [ m are the angular eigenvalues and eigenfunctions which in 
general can only be obtained perturbatively. At this point, the perturbation of the Kerr 
black hole has become a typical Quantum Mechanics work in which one has to solve a 
Sturm-Liouville problem for two coupled radial and angular ODE's. The Teukolsky Master 
Equation can be cast in a more compact form by introducing a "connection vector" whose 
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components are: 



j-ir 

r<t> 



1 



M{r 2 - a 2 ) 



A 



[r + ia cos ( 



— (r-M) 



It's easy to prove that: 



1 



1 



a(r — M) cos # 
-+^- 



A 



1 



r^r„, = — cot 2 e + m£ 



and consequently the Teukolsky Master Equation assumes the form: 

[(V" + sr")(V M + sFA - As 2 ^]^ s) = AtxT 



(4-12) 



(4-13) 



(4-14) 



where ^ is the only non vanishing NP component of the Weyl tensor in the Kerr background 
in the Kinnersley tetrad (4-4). Equation ( 4-14)) gives a common structure for these massless 
fields in the Kerr background variating the "s" index. In fact, the first part in the lhs 
represents (formally) a D'Alembertian, corrected by taking into account the spin-weight and 
the second one is a curvature (Weyl) term linked to the "s" index too. This particular form 
of the Teukolsky Master Equation forces us to extend this analysis in the next sections. 



§5. Generalized wave equations 

The compact form of Teukolsky Master Equation we have shown suggests a connec- 
tion between the perturbation theory and a sort of generalized wave equations which differ 
from the standard ones by curvature terms. In fact generalized wave operators are know 
as De Rham-Lichnerowicz Laplacians and the curvature terms which make them different 
from the ordinary ones are given by the Weitzenbock formulas. Mostly known examples in 
electromagnetism are 

• the wave equation for the vector potential A^©: 

V a V a A ll -R*A x = -4xJr , V a A a = (5-1) 

• the wave equation for the Maxwell tensor 0*: 

VVj.PU + R p ^xF^ - R p \F vp + R p v F Xp = -8nV [p J u] (5-2) 
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while for the gravitational case one has 

• the wave equation for the metric perturbations (gravitational waves)©: 

V a V a h^ v + 2R a / 3flI ,h al3 — 2R a ^h u ) a = 0, 

Vji" = 0, = hp, - -g^h a a (5-3) 

• the wave equation for the Riemann Tensor© (in vacuum for simplicity) 

V M V M -R a/ 3 7 5 — Ra/3p(rR^S f " 7 ~ 2(R apya Rf3 P f — RapSaR^-y) = . (5-4) 

These equations are "non minimal," in the sense that they cannot be recovered by a minimal 
substitution from their flat space counterparts. A similar situation holds in the standard 
Quantum Field Theory for the electromagnetic Dirac equation. In fact, applying for instance 
to the Dirac equation an "ad hoc" first order differential operator one gets the second order 
Dirac equation©'©' 



(ifi — eA + m) (ip — eA — m)ip 



e 



{id ll -eA lM ){id»-eA»)-~a» v F l 



pv 



m 2 



V = 0, (5-5) 



where the notation is obvious. It is easy to recognize in eq. fl5-5|) a generalized Laplacian and 
a curvature (Maxwell) term applied to the spinor. Moreover this equation is "non minimal" , 
in the sense that the curvature (Maxwell) term cannot be recoverered by electromagnetic 
minimal substitution in the standard Klein-Gordon equation for the spinor components. 
The analogous second order Dirac equation in presence of a gravitational field also has a non 
minimal curvature term and reduces to the form: 

(V a V Q + m 2 + ^R)iP = (5-6) 

Quite recently it was shown by Mohanty and PrasannaH), that the Maxwell tensor wave 
equation ( |5-2| ) could be framed into the context of the QED vacuum polarization in a back- 
ground gravitational field. They start from a classical work of Drummond and Hathrellil) 
which derives higher derivative coupling arising from the QED loop corrections to the 
graviton-photon vertex, encoded in the Lagrangian, valid for u 2 < a/m 2 (a is the fine 
structure constant and m e is the electron mass): 

C = --F^ + aRF^ + bR^F\ + cR^F^F" (5-7) 

with the coefficients a = — ^g-, b — c = — ^g-- The physics of these "tidal" terms 

is that the effect of virtual electron loops gives to the photon a typical size proportional to 
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the Compton wavelenght of the electron A c . The field equations obtained from the modified 
lagrangian fl5-7|) are: 

V^F^ + (T = (5-8) 

where Q v represents the non linear quantum correction the explicit (long) form of which is 
given ini3 ) . Manipulating the previous relation one gets the wave equation: 

+ {R^xF"" + R\F vp - R p u F Xp } + [V U (2 X - V A ^] = . (5-9) 

In the range of validity of the effective action the Riemann and Ricci terms in the curl 
brackets of ( |5-9| ) which are present in Einstein's gravity are larger than the terms in the 
square brackets, arising from the loop corrections, and consequently it is clear that the 
classical second order wave equation for the Maxwell tensor describes with good accuracy 
this quantum phenomenon too. 



§6. Vacuum gravitational fields 

The Weyl and the Riemann tensor are connected by the relation: 

-~(g ai g(3s - g a sg/3y)R , (6-1) 

In a vacuum spacetime, R a p = 0, R = and consequently the Riemann tensor coincides with 
the Weyl one. The Weyl tensor in vacuum satisfies the generalized de Rham wave equation 
(the vertical bar indicates the covariant derivative) 

^(dR)Cai 020304 ^01020304 |6 (C a\a2^bcazai ~\~ 4C ai [GyCaslcb,^) . 

(6-2) 

The name "generalized de Rham wave equation" comes from the nice compact form which 
can be given to this equation 

AdR)^! a 2 a 3 a 4 = (SD + D5)C ai a 2a , iai = , (6'3) 

by using the divergence operator 5 and the covariant exterior derivative D^S'Ei} (the 
divergence operator 5 and the covariant exterior derivative D must not be confused with the 
NP frame derivatives denoted by the same symbols). This is the most natural generalization 
to tensor-valued forms of the ordinary de Rham wave operator 

A m = 6d+d6, (6-4) 
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whose action is defined only on ordinary differential forms. A detailed study of the general- 
ized de Rham operator will be discussed elsewhere 0. 

The explicit form taken by [Z^(dR) C] a b c d = in the Newman-Penrose formalism is now 
derived and simplified somewhat using the commutator identities to eliminate some second 
derivatives and the Bianchi and Ricci identities to eliminate some coupling terms between 
different Weyl tensor scalars. 



[Z\(dR)C]i3i3 = (wave equation for \Pq) 



Recalling that $b = _ C1313, select the 1313 component of eq. (J6 

[-DA -AD + 5*5 + 55* + (9 7 + 7* - y - y*)D + (7e - e* + p* + p)A 
+ ((3* + n - 9a - t*)5 + (n* - t - 7/3 - a*)5* + 4(£> 7 + Ae - 5a - 5*f3) 
+ A(aa* - p*7 + r* (3 + tic - (3* (3 - 87c - e 7 * + 7 e* - an* + ye - (3n - a\ 

(6-5) 

+ kv + ra + 8a/3 - py - p 7 + y*e + 3& 2 )}% + 4[2p5 + 2a5* - 2tD - 2k A 
+ 5p — Dt + 5*a — Ak + 71* p — pa* + 5re — re* — 5pj3 + rp* + an — Ky* 
- Ky - 7aa + k 7 * + 7k 7 - ar* + a(3* - W^/ x + [2A(ap - kt)}& 2 = . 

This equation only contains second order directional derivatives of which can be rewritten 
first using the identities 



DA + AD = 2D A + [A, D] 

55* + 5*5 = 255* + [5*, 5] (6-6) 



and then the commutators in these expressions can be replaced by first order derivatives 
using the commutator identities fl2-7|) i and ( p-7D 4 , in order to reduce the number of second 



order differential operators acting on 
Equation (6-5) then becomes 

[-2DA + 255* + (87 - 2y)D + (6e - 2e* + 2p*)A + (2vr - 8a)5 
+ (2tt* - Qf3 - 2a*)5* + 4(£> 7 + Ae - 5a - 5*f3) + 4(aa* - p* 7 
+ t*I3 + ttt - (3*3 - 876 - e7* + 7e* - an* + ye - I3n - aX 

(6-7) 

+ kv + ra + 8a(3 - py - p>y + y*e + W 2 )}% + 4[2p5 + 2o5* - 2tD - 2kA 
+ 5p — Dt + 5* a — Ak + n* p — pa* + 5re — re* — 5p/3 + rp* + an — Ky* 
- Ky - 7aa + k^* + 7kj - ar* + a/3* - + [24(crp - kt)\P 2 = 

Next the Bianchi identities can be used to simplify this expression and introduce terms 
involving \P to decouple it whenever possible. This task is achieved by the elimination of 
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two directional derivatives of #1 in eq. (6-7) using two Bianchi identities 

#i3[i3|4] = : D&! - 5*% + (4a - tt)% - 2(2p + e)^i + 3k& 2 = , 
i2i3[i3|2] = : 6$! - A% + (4 7 - n)V - 2(2r + (3)^ + 3a& 2 = . (6-8) 

Equation (6-7) then becomes 

[-2£>Z\ + 255* + 2(4 7 - p)D + 2(3e - e* + p* + 4p)Z\ + 2(tt - 4a)<5 
+ 2(tt* - 3(3 - a* - 4r)<T + 4(£> 7 + Z\e - 5a - 5*13) + 4(aa* - p* 7 
+ r*/3 - T7T - tf*/? - 8 7 e - e 7 * + 7 e* - avr* + pe - (3n - aX 

(6-9) 

+ kv + 9ra + 8a/9 + pp - 9p 7 + p*e + W 2 )]% + 4[2cr5* - 2/tZ\ 

+ 5p — Dt + 5*o" — Ak + 7r*p — pa* + re — re* — p/3 + rp* + o"7r — /tp* 

- Kfx- 7cra + /t 7 * + 7/t 7 - or* + a"/?* - 3«Pi]«?i = . 

Finally this (exact) equation can be given a Teukolsky-like form which is very helpful in 
type D geometries where most of the known results have been derived. To accomplish this 
additional step, multiply the equation by —1/2 and use three Ricci identities in two different 
ways: a) in (6-9) first use the identities 

[l/2(#i2i 2 - #3412)] : At = D 1 - a(r + tt*) - (3{r* + vr) + 

+ 7 (e + e*) + e( 7 + 7*) - rvr + uk - & 2 , 

6-10 

[1/2(^1234 - #3434)] : S*(3 = 5a- (pp - \a) - aa* - (3(3*+ 

+ 2a(3 - 7 (p - p*) - e(p - p*) + V 2 ; 
b) rewrite the other Ricci identity [#2431] in the following form 
[#2431] : Dp-5n- (p*p + o-A) - 7r(7r* - a* + (3) + p(e + e*) + uk - & 2 = , (6-11) 

multiply it by \P and c) add this (which is actually zero) to the wave equation. The result 
is 

[(D - 3e + e* - 4p - p*)(A + p - 4 7 ) - (5 + tt* - a* - 3(3 - 4r) 
(5* + n- 4a) - 3W 2 }% + 3(<t A - kv)% + 2[2kA - 2o5* - 5*a + An 

(6-12) 

— 5p + Dt — 7/c 7 — (3* a + 7aa + r*cr + up* + up — pix* + re* + (3p 

- er - k 7 * -1x0 - rp* + pa* + Wi\^i = ; 
d) finally use in eq. (6-12) two other Ricci identities 

[#3143] : Sp - 5* a = p(a* + (3)- a(3a - (3*) + r(p - p*) + n(p - p*) - ^ , 
[#1312] : Dt-Ak = p(r + tt*) + a(r* + tt) + r(e - e*) - «(3 7 + 7*) + 3^. (6-13) 
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The final form for eq. (6-12) is then 

[(D - 3e + e* - 4p - p*)(A + fi - 4 7 ) - (5 + tt* - a* - 3/3 - At) 

(5* + 7i- Aa) - 3^ 2 ]^o + 3(o- A - kv)W q + [AkA - 4o<f - 45V + AAk 

- 20/t7 - 4/9*0- + 20o-a + 4rV + 4k//* - 4/t 7 * + 10#i] #i = , 



(6-14) 



which is the exact form for the same equation for \Pq found by Teukolsky in the very special 
case of a type D geometry in the first order perturbation regime (see below). Moreover 
this equation coincides with that obtained by Stewart and Walker© following a different 
approach and in the Geroch-Held-Penrose (GHP) formalism© 



( p'p - d'd - pp - 5pp' + fd + 5rd' + Aaa' - Akk' - lOS^o 
+ [4p'« - 4<9'o - 4(p' - 2p')n + 4(f - 2r')o + 10#i]#i 
+ (-4op + And - YIkt + I2pa)$ 2 = . 



(6-15) 



The equivalence of (6-14) and (|6 • 1 5|) is now explicitly shown. Recall the action of the 
two-index {p, g}-operators p, p', d, & on two-index {p, g}-objects. Taking into account the 
weights of the ^„'s (i.e., {p, q} = {A — 2n, 0}), those of the spin coefficients and the definition 
of the weighted derivative operator of an object r\ (which increases its {p, q] indices) 



p?7 = (D — pe — qe*)rj, 
p'p = (A-jry - qj*)rj, 
Or] = (5 — pj3 — qa*)r], 



increasing of type { + 1 , + 1 } , 
increasing of type { — 1,-1}, 
increasing of type {+1,-1}, 
increasing of type { — 1 , + 1 } , 



(6-16) 



d'i] = (5* — pa — q/3*)r], 
eq. ( |6-15|) becomes: 

[(A - 57 - j*)(D - 4e) - (5* - 5a + (3*)(5 - 4/3) + p*(D - 4e) 

- 5p(A - 4 7 ) + r*(5 - A(3) + 5t(S* - Aa) - A(a\ - ku) - lOS^o 
+ (-AaD + Ak5 - \2kt + \2pa)V 2 + [A(A - 37 - j*)k 

- A(S* -3a + (3*)a - A(-p* + 2p)K + 4(r* + 2n)a + lO^i)]!^ = . 



(6-17) 



where the second order operator has to be rewritten via the commutation rules (2-7 ). 

Then, by using the two Ricci identities l/2[i? 1212 — i?34i 2 ] and l/2[i? 1234 — -R3434] (see 
eq. (6-10)), adding equation [i? 2 43i] multiplied by \P and using the two Bianchi identities 
#i3[43|2] = and #i3[2i|4] = 0, namely: 



A$x - 5^2 - v&o - 2(7 - (J,)&i + 3r^ 2 - 2o^ 3 = , #i 3 [43|2] = 
5*& x - D&2 - A^o + 2(tt - a)$ x + 3p& 2 - 2k& 3 = , #i 3 [2i|4] = 



(6-18) 
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one gets exactly eq. (6-14). 



^(dR)C]2424 — 



(wave equation for #4) 
The equation for ^4 can be immediately obtained from the equation for & using the ex- 
change symmetry 1 <-> n, m <-> m* of the Newman- Penrose quantities. Under this operation, 
eq. (6-14) becomes: 

[(A + 3 7 - 7* + 4p + fi*)(D + 4e - p) - (8* - r* + (3* + 3a + An) 

(8 - r + 4/3) - 3^ 2 ]^ 4 + 3(aA - ki/)#4 + [4A5 - 4zaD + 45A - ADv (6-19) 

- 20z/e - 4a*A + 20A/3 + 4vr*A + Avp* - Ave* + 10^ 3 ]^ 3 = . 

Working in the GHP formalism, the corresponding relation found by Stewart and Walker 
0) is: 



[( p'p - a'a - (V + p')p - Pp' + ( 4t ' + ^)<3 + r<9' + 4pp' - 4rr' - 2# 2 ]# 4 

+ [4p«' - Ada' - 4(p - 2p)«' + 4(f - 2r)o / + 10^ 3 ]^ 3 (6-20) 

+ [-4</p' + Ak'8' - Uk't' + 12pV]# 2 = 0. 



[A( d R)C 



1213 



(wave equation for W\) 



Selecting the component C1213 ($1 = — C1213) in eq. ( |6-2|) one has 



[88* + 8*5 - DA - AD + (7* + 5 7 - p - p*)L> + (p + p* + 3e - e*)Z\ 
+ (tt + (3* - t* - 5a)8 + (n* - 3(3 - r - a*)8* + 2(L>7 + At - 8a - 8* '(3 
+ 5r7r + aa* — 47e — /ry — p*7 + p*e — 67* — (3(3* + -ye* — /?7r 

- 5(tA + /3r* - an* + pe + 4a/3 - 5pp + ra + 5ki/ - 3& 2 )}&i 

+ [2uD + 2nA - 2X8 - 2p8* + Du + An-8X- 8*p + r*p - up* + rX (6-21) 
+ vr/i* + Aa* - 777T - 5z/e + 7ap - tt7* + z/e* - z/p + 5 X(3 - Xn* - p/3* + Q& 3 }% 
+ 3[-2tL> - 2kA + 2p<5 + 2(T<5* - Dr - An + 8p + 5V + /3V - /3p - ar* 

— Kpf + n-f* — 3cra — re* + er — a*p + rp* + 37K — p/t + an + pir*]W 2 
+ 12(ap - kt)& 3 = . 

Using the commutation rules ( |2-7| ), multiplying by —1/2, and using the six Bianchi 

identities i£l3[13|4] = 0, #13[13|2] = 0, -Ri3[21|4] = 0, -Rl 3 [43|2] = 0, #42[13|2] = i?42[13|4] = 

and the five Ricci identities [#3143], [-R1312], [-R2443], [R2421], [5CR1234 - Ram))], eq. (6-21) 
becomes: 
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[DA - 55* + 2(/i - 7 )L> + (-3p - p* - e + e*)A + 2(a - n)5 + (p + a* + 3r - n*)5* 

— I?7 — Z\e + 25a + 2p*7 — 7e* + «7r* + kz/ + e 7 * + 6/ry + 4 7 e — 6pp — 7ra + 5r7r — [3t* 

- 2aa* - 2a(3 + 3pir - 4pe + 7& 2 ]&i + [A<J - ^ + 5\ - Dv - 3rA +ue- a*X (6-22) 
+ 3pv - vt* + Att* + up* - f3\ - 2$> 3 }% + 3[k5 - oD + 2<xp + 2k/3 - 2kt - 2(xe]^ 3 

+ 3[Z\k - 5* a - 2fiK + 2an - af3* + Kji* - 3k 7 + ar* - k 7 * + 3aer]#2 = 0. 



Now adding to the previous relation the equation [-R2431] (see eq. (|6-11|) ) multiplied 2$i and 
using another Ricci identity, [|(i?i2i2 — -R3412))], (6-22) takes the compact form: 



[(D — e + e* — 3p — p*)(A + 2p - 2 7 ) - (5 - a* - (3 + n* - 3t)(5* + 2tt - 2a) 
+ 6^2]$! + 2(«i/ - <xA)#i + [A5 - vD + 5\- Dv - 3rA + i/e - ct*A + 3pz/ 
- ue* + Att* + up* - (3X - 2& 3 }% + 3[k5 - oD + 2ap + 2k/3 - 2kt - 2ae}& 3 
+ 3[Z\k - 5*a - 2fiK + 2air - af3* + Kp* - 3k 7 + ar* - k 7 * + 3aa}& 2 = . 
The analogous equation in the GHP formulation is 

[(p - p - 3p)(p' - 2p') -(d-f 1 - 3r){d' - 2t') + 2{aa' - ««') + 6^]^ 
+ [3k& - 3<xp + 6a p - 6kt}& 3 

+ [p«/ + k'\> - a'd - da' + a'f' + 3aV - Kp - 3k p - 2& 3 }W 
+ [3pK - 35V + 6kp' - 6ar' + 3af - 3«p']^ 2 = . 



(6-23) 



(6-24) 



^(dR)C]l242 — 



(wave equation for \P 3 ) 
\P 3 can be easily obtained, as was done for 1^4, applying the 



The equation for C1242 = 
interchange symmetry 1 <-> n, m <-> m* to (6-23), yielding: 

[(Z\ + 7 - 7* + 3p + p*){D -2p + 2e) - (5* + /3* + a - r* + 3tt)(<5 - 2r + 2(3) 
+ 6^ 2 ]^3 + 2(ki/ - cxA)^ + [kZ\ - cxcT + Z\k - 5V - 3cttt + ^7 - /3V + 3pK 
- kj* + or* + up* -aa- 2^i]^ 4 + 3[\A - u5* + 2A/i + 2va - 2un - 2\j)&\ 
+ 3[5\ -Dv- 2pv + 2Ar - \a* + up* - 3ue + Att* - vt* + 3/3A]^ 2 = . 
The analogous equation in the GHP formulation is 



(6-25) 



W -p'- 3p')(p - 2p) -(d'-f- 3r')(d - 2t) + 2{a'a - k'k) + 6^3 
+ [3k'& - 3a'p' + 6a' p' - Qk't'}^ 

+ [Jj'k + 4,' - 0-5' - 5V + o-f + 3ar' - np' - 3np' - 2& 1 }& 4 
+ [3pK' - 35a' + 6k' p - 6a V + 3ar ' - 3k'p)& 2 = . 



(6-26) 
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[^(dR)C]l342 — 



(wave equation for ^2) 
Selecting the component 1342 of (|6-2|) ($2 = -^1342) one has 

[-DA -AD + 55* + 5*5 + (7 + 7* - p - p,*)D + (p + p* - e - e*)A 

+ (n + (3* -a- r*)5 + (vr* + (3 - a* - t)5* + 12ttt - 12pp + \2vk - 12Xa - 6^2 

+ \AvD + An A - AX5 - Ap5* + 2Dv + 2 An - 25X - 25* p + 2pr* - 6^n + 2ue* - 2p*v 

- 2p(3* + 2A/3 + 2\a* + 2np* - 2^* - 2ue + 2rA - 2pv - 2\n* + Gap + 4^]^ (6-27) 
+ [-ArD - AkA + Ap5 + 4o-5* - 2Dr - 2Ak + 25 p + 25* a + 2n*p - 2e*r + 2rp* 

- 2pn + 2aa + Qp(3 - 2t*o + 2a(3* + 2ttct - 2np* - 2k-j - 6er - 2pa* + 2k7*]^ 3 
+ [Apa - 4tk]^4 + [A\p - Avtx + 2$^% = . 

Using the commutation rules ( |2-7| ), multiplying (6-27) by —1/2, and using the four Bianchi 
identities i?i 3 [ 2 i|4] = 0, -Ri 3 [43|2] = 0, -R 42 [i3|2] = -R 42 [i3|4] = 0, the four Ricci identities [-R3143], 
[-R1312], [-R2443], [-R2421] and adding eq. ( |6-11| ) once multiplied by 3$2, relation (6-27) takes 
the form: 



[(D + e + e* - 2p - p*)(A + 3p) - (5 + vr* - 2r + (5 - a*)(5* + 3tt)}W 2 
+ (3A<r - 3z/k)^ 2 + 2\-uD + X5- Du + 5\ + up* + \n* - 2r\ 

- Xa* + X(3- ve* -P6+ 2pv + ^3]^ + 2[kA - a5* + An 

- 5* a + aa + Kff + 2pn + r*a - kj* -jk- 2na - af3*]& 3 - $b$4 = 

The corresponding equation in the GHP formalism is 



(6-28) 



[(f> -2p- p)$ - 3p') - (d - 2t - f')(& - 3r') + 3(kk' - aa')]^ 2 

+2[«/fc + \>k' - a'd - da' + 2ra' + a'f' - 2 P k' - k'p]$i (6-29) 
+2[k$>' + b'« - d'a - ad' + ^1 + 2rV - 2p 'k + f a - np'}^ 3 - ^ ^4 = . 

At this point the complete set of (exact) wave equations satisfied by the Weyl scalars has 
been given. We can now perform their perturbation expansion. 

6.1. The perturbation theory 
Introducing the operators: 

t = [(D - 3e + e* - Ap - p*)(A + p - 4 7 ) - (5 + vr* - a* - 3(3 - At)(5* + vr - Aa) - W 2 ], 
W = 3(aA- ku), (6-30) 
Vi = [AkA - Aa5* - A5*a + AAk - 20/«7 - 4/3 V + 20act + At* a + Anp* - 4/t 7 * + 10#i], 
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where the letter % here stands for the operator on $o found by Teukolsky in the perturbation 
regime; this is a notation extending that due to Lousto and Campanelliil'. In the same way 
a generic operator x acting on will be denoted by Xi- Thus, eq. (6-14) can be written as: 

(f + W )% + V&x = . (6-31) 

Performing the standard perturbation expansion of the various quantities a la Teukolsky: 
% = # (°) + ^ Q (1) + + a = <rC°) + a® + ...,D = D<® + £>« + ... and so on, one finds 

[(T (0) + f (1) + f (2) + ...) + (Hf + + W {2) + -W^ + ^o (1) + + 7 ) 
+ {V? + V? + V? + ...)(^ (0) + + *i 2) + •••) = , 

giving an inductive procedure for obtaining a hierarchy of equations of the following type: 
(f {0) + W (0 Vo (0) + VM 0) = 0, 

(f (1) + W «)tff + (t (0) + VpfW 3 + (A (0) ^i (1) + A (1) ^i (0) ) = 0, 



(6-33) 



(f {n) + n^Vf + (r (0) + vvf y (n) + {vtM 0) + vfM n) ) = s [ o ] , 



where 



n-1 

= - Y, [(% (n ~ P) + yvt~ P) )^o P) + Vt~ p) v[ p) ] , (n > 2) . (6-34) 
P =i 

For every n > 2 perturbation equation, the "source term" Sq^ is composed of quantities at 
most of order (n — 1), which in principle are known. 

For a type D background geometry it is well known that 



K C°) = £r(°) = A<°> = = 0, 



(6-35) 



so VVq ^ = Wq 1 ^ = V[ QS) = and the recurrence scheme (6-33) becomes: 



= 0, 

r (0) y/ (1) 



-r(0),7v(n) _ c-(n) 



(6-36) 
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One can study the perturbation of eq. (6-19) exactly in the same way as before. In fact, 
defining the operators 



% = [(A + 3 7 - 7* + 4/i + p*){D + 4e - p) - (5* - r* + /F + 3a + 4vr) 
(8 - t + 40) - 3V 2 ], 

~ (6-37) 
V 3 = [4X5 - 4vD + 45 A - 4Dv - 20ue - 4a* A + 20A/3 + 4vr*A + 4up* - 4ue* + 10# 3 ], 

>V 4 = 3(crA - ku), 

(where the letter 71 here stands for the operator on $4 found by Teukolsky), eq. (6-19) can 
be written as 

(t 4 + W 4 )^4 + V 3 & 3 = . (6-38) 
The hierarchy of equations in this case is: 

(fi 0) + wi 0) )^ 0) + pf^ 0) = o, 

(t 4 (i) + wivf + (t 4 (o) + wi°y 4 (i) + (pf^ 3 (1) + = 0, 

(6-39) 

(f 4 (n) + vPi n Vi 0) + & 0) + vPi Vi n) + (Pfty 0) + ^ 0) ^ (n) ) = $i n \ 

where 

n— 1 

5| n) = - £ [(T 4 (n ~ p) + Wi n ~ p V 4 (p) + pf^Vf] , (n > 2) . (6-40) 
P =i 

In the case of a type D geometry w|°' ) = = = and the set of equations (6-39) 

becomes: 



= 0, 



(6-41) 



The previous results are an alternative but equivalent formulation of the work of Lousto and 
Campanelli concerning the higher order perturbations of a type D vacuum metric. The same 
perturbation formulation can easily be written for eqs. (6-23), (6-25) and (6-28) allowing a 
systematic formulation of the gravitational perturbation for every vacuum spacetime. In 
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fact by introducing the operators: 

Ax = [(D — e + e* — 3p — p*)(A + 2/i - 2 7 ) - (5 - a* - (3 + n* - 3r)(5* + 2rr - 2a) + 6& 2 ), 
B x = 2(ku-(jX), 

K = [\S - U D + 5X-Du- 3rA + ue - a*\ + 3pu - ue* + Xtt* + up* - f3X - 2# 3 ], (6-42) 
Q 3 = 3[k5 -aD + 2ap + 2k/3 - 2kt - 2ae], 

M 2 = 3[Ak - 5* a - 2pK + 2an - a/3* + Kpf - 3^7 + ot* - K/y* + 3aa], 
eq. (6-23) can be rewritten as: 

(Ax + Bx)#x + iC % + Qs^s + -M 2 ^ 2 = , (643) 

and one can develop a hierarchy of perturbation equations, analogously to what we did 
above. In the case of the first order perturbation of a type D background one finds 

jgjW = BP = = Qf ] = Mf = 0, (644) 

so that the linearization of eq. ( 643|) gives 



A { M 1] + M { M 0) = 0, (645) 

which unfortunately couples Weyl scalars and perturbed spin coefficents. Relations for W% 
can be worked out in exactly the same way. For #2 the following operators are introduced: 

T 2 = \(D + e + e* - 2p - p*)(A + 3p) - (5 + n* - 2r + (3 - a*)(5* + 3rr)], 
Q 2 = (3Acr - 3un), 

_ (646) 
M = 2[-zaD + X5 - Du + 5A + up* + Att* - 2rA - Xa* + X(3 - ue* - ue + 2pu + $ 3 ], 

C 3 = 2[kA — a5* + Ak — 5* a + aa + np* + 2p,K + t*o — kj* — — 2-rij — a ■/?*]; 
eq. (6-28) can then be rewritten as: 

(T-i + £2)^2 + M^i + £3^3 - « = . (647) 
Again the perturbation theory can be directly developed. In particular, for the type D case 

Qf = §P = A/f = 4 0) = , (648) 
and the first order perturbations become: 

+ (;F 2 (1) + gP)^ 0) = . (649) 
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Thus, even in the very special case of a type D geometry, the perturbation equations for 
$i,$2 and $3 are coupled — order by order — with other perturbation quantities. This is due 
to the presence, in the exact equations, of terms involving ^ which are responsible for the 
appearance of the perturbations of the spin coefficents and of the directional derivatives. 
To obtain decoupled equations for these three Weyl scalars in the perturbation regime, in 
principle, one has to define linear operators with derivatives of order higher than the second. 

We conclude this section by pointing out that our perturbation relations (6-36) and (6-41) 
assume a more compact form if compared with those found by Lousto and Campanelli be- 
cause of the completely different approach and the number of simplifying operations (always 
explicitly specified) that we perform at the level of the exact theory. They use 'background 
operators' acting on the exact Bianchi identities and then linearize the results order by order 
(see also their remark after eq. (6) in©), following the Teukolsky prescription. The way of 
using exact 'ad hoc' operators acting on the exact Bianchi identities and then linearizing the 
results was introduced by Stewart and Walker but only for first order perturbations. 

By extending such a prescription to any order one obtains our relations. The de Rham 
Laplacian mediates this task at the level of the exact theory, also explaining the nature (i.e. 
exact 'wave equations in curved spacetime' for some of the NP Weyl and Maxwell scalars) 
of the relations found algebraically by Stewart and Walker. 



§7. Vacuum spacetimes: test electromagnetic fields 

In this section we study the electromagnetic field equation in the absence of sources, 
J M = 0. This choice has to be related to the "exact" (in the sense of the field equations) 
nature of our approach, which goes beyond the perturbation level and would be corrupted 
by the introduction of a source term without the automatic appearance of a corresponding 
term in the Riemann tensor. On the other hand, eq. ( |7-1| ) below is itself an approximation, 
because the exact field equations should include the Riemann and Ricci terms created by the 
EM field and should be coupled to the gravitational de Rham Laplacian. The formulation 
presented here is based on the hypothesis of test fields (no backreaction) and the goal of 
this analysis is to recover the existing perturbation theory, including the work of Teukolsky, 
Fackerell and Ipser and many other authors. 

The sourcefree Maxwell equations 5F = = dF for the Maxwell 2-form F (representable 
as 4 complex equations for the associated NP scalars <po, <f>i, 02, namely eqs. (330)-(333) of 
Chapter 1 o 

fi) 

) immediately imply the vanishing of the de Rham Laplacian of the Maxwell 
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2-form. In a vacuum spacetime, this takes the form 

[A {dR) F] ab = W c + c ahcd F cd = o 



(7-1) 



Expressing these equations in the Newman- Penrose formalism, one has the following wave 
equations. 

( wave equation for O ) 



13 







Studying the component "13" of ( [7T| ) we get the equation for O : 

[DA + AD- 5*5 - 55* + (-5 7 - 7* + p + fi*)D + (-3e + e* - p* - p)Z\ 
+ - vr + 5a + t*)5 + (-7T* + r + 3/5 + a*)5* + 2(5*/? + 5a - £7 - Z\e 
+ 7*6 - 7e* - r7r + 7p* + n*a + A^e - (3r* + pp + (3(3* - ep* - aa* - Aa(3 
+ p7 - kv + /far + aX - ra - fie - ^ 2 ]0o + [4(rL> + kA - p5 - a5*) 
+ 2(Dr + Ak- 5p- 5*a + e*r - pn* - er + err* - ^7* - rp* - a(3* + pa* 
+ p/t + p*/t + /3p - on - 3^7 + 3aa + 2#i)]0 1 + [4kt - 4ap - 2<Z' O ]02 = . 



(72) 



First we use the commutation rules ( |2-7| ); the next step is to divide the resulting equation 
by 2 and to use the two Maxwell equations: 



(D - 2p)(f) 1 - (5* + 71 - 2a)0 o + ^2 = 0, 
{5 - 2r)0! -(A + n- 2 7 )0 O + vfc = 0, 



(7-3) 
(74) 



to eliminate two of the directional derivatives of <\>\. Then using four Ricci identities: 
[l/2(i?i2i2 - i?34i 2 )] : Ae = D 1 - a(r + vr*) - (3(r* + vr) + 7 (e + e*) 

+ e( 7 + 7 *) T7T -|- VK <^ 2) 

[l/2(-Ri234 - ^3434)] : 5*/? = 5a - (pp - \a) - aa* - (3(3* + 2a(3 

- 1(P ~ P*) - <H - A 4 *) + #2, 



(75) 



together with the two in (6-13), rewriting the Ricci identity [-R2431] as in ( |6-11| ), multiplying 
it by 0o and adding the resulting quantity (which is actually zero) to the wave equation one 
obtains: 

[(D — e + e* — 2p — p*)(A + p - 2j) - (5 - (3 - a* - 2r + tt*)(5* + tt - 2a)]0 o 
+ (a A - Kv)<f> + 2[kA - (x5* + zAk - 5V - k 7 * - a/5* + Kp* + cxr* + 3aa (7-6) 

- 3^7 + 2^i] 01 -^002 = 0, 

where in the first line it is easy to recognize the operator found by Teukolsky for the elec- 
tromagnetic perturbations of a type D metric. One can also verify that (7-6) is the electro- 
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(7-7) 



magnetic equation of Stewart and Walker in the GHP formalism: 

[pp' - d& - p'p - (2p + p)p' + r'd + (2r + f>)& + 2pp' - 2rr' + ^ 2 ]0o 
+ [2p're - 2d' a - 2k{p' - 2p') + 2a{f - 2r') + 4lPi]0i + [2/t<9 - 2<xp 
+ 2pa - 2kt - ^o]0 2 = . 

In fact "translating" (7-7) into the standard Newman-Penrose form one gets: 

[{D - e + e*)(Z\ -2y)-(6-l3- a*){5* - 2a) + p(D - 2e) 
- (2p + p*)(A - 2 7 ) - tt(5 - 2/3) + (2r - tt*)(5* - 2a) - 2pp + 2rvr + ^ 2 ]0o 
+ 2[(A - 37 - 7*)k - (5* -3a + (3*)a + ft(p* - 2p) + <x(r* + 2tt) + 2^]^ 
+ [2k(5 + 2(3) - 2a{D + 2e) + 2pa - 2kt - %}<p 2 = . 

Using the following two Maxwell equations: 

{D-p + 2e)0 2 - (5* + 2tt)0i + A0 O = 0, 
(5-t + 2/3)0 2 - (Z\ + 2p)0 x + z/(/) = 0, 

in the previous equation and adding the Ricci identity [-R2431] (see eq. |6-11|) multiplied by 
0o, one finds identically eq. (7-6). 

(wave equation for 2y ) 



(7- 



(7-9) 



A(<m)F) 2 4 = 



The equation for 2 can be immediately obtained form eq. (7-6). Using the exchange 
symmetry 1 <-> n, m «-> m* of the Newman-Penrose quantities and the relation (|2-18| ) 

gets 



one 



- {[(A + 7 - 7* + 2/i + //)(£> - p + 2e) - (<f + a + (3* + 2tt - r*)(5 - r + 2/?)]0 2 
+ (o-A - Kv)(f) 2 + 2[A<5 - vD + 5X - Dp - ue* - Xa* + up* + Att* + 3/3A (7-10) 

- 3z/e + 2^ 3 ]0i -^40o } = . 

The overall minus factor in our case can be eliminated but this form is convenient for further 
developments in the presence of sources. The GHP analogue of this equation is 



[p'p - d'd - pp' - (2p' + p') P + r& + (2r' + f)d + 2p'p - 2t't + <^]02 

+ [2p«/ - 2da' - 2k' (p - 2p) + 2a' (t' - 2r) + 4<Z' 3 ]0 1 + [2k'& - 2a'p' (7-11) 

+ 2p'a' - 2k' t' - # 4 ] O = 0. 



[AdR)^]i2 - [A (dR) F] u = 



[wave equation for <pi) 
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To obtain the equation for 0i one has to subtract the components "12" and "34" of (|7-1|) 
in order to eliminate all the 0* terms. The resulting wave equation, divided by a factor of 
2, is: 

[DA + AD- 5*5 - 55* + (-7 - 7* + fi + fi*)D + (e + e* - p* - p)A 

+ {-j3* - n + a + t*)5 + (-7T* + t - (3 + a*)5* + 4(pp + Act - rvr - vk + ^ 2 )]0i 

[-2zaD - 2nA + 2X5 + 2u5* - Dv - An + 5X + 5*u - \a* + pv + 7r 7 * - ur* 

(7-12) 

- tA - 3a/i + Att* - 717/ + - A/3 - z/e* + z/p* + 3771 + ve - 2<Z' 3 ]</> 

+ [2tD + 2/tZ\ - 2p5 - 2a 5* + Dt + Ak- 5p- 5* a -an- rp* + 7 k - k 7 * 

+ 3re + tV - 3p(3 + p*K - aa + re* - a/3* + pa* + pLK- pn* - 2^i}<p 2 = 0. 

The procedure to simplify the previous relations is the usual one. The starting point is 
to rewrite the second order part by using (6-6); dividing the resulting equation by 2 and 
using all four NP Maxwell equations, some of the directional derivatives of 0o an d 02 are 
eliminated in favor of the derivatives of 0i. Now, using the Ricci identies: (6-10), (6-13) and 

[#2443] : SX - 5* pl = v(p - p*) + - pl*) + n(a + (3*) + A(a* - 3(3) - ^ 
[Rmi\ ■ Dv- An = fi(n + r*) + X(n* + r) + tt( 7 - 7*) - v(3e + e*) + ^3 ; 
adding the Ricci identity fl6Tl|) multiplied by 20i, we finally obtain 

[(D + e + e* - p - p*){A + 2p) - (5 + (3 - a* - r + n*)(5* + 2vr)]</> 1 

+ [-vD + X5 - Dv + 5X- Xa* - e*v + pv - r X + A/3 - ve + vp* + tt*A]0 o (744) 

+ [kA — a 5* + Ak — 5* a — a (3* — K7* + pin — an + up,* + ar* + ao — K7]02 — . 

It is easily seen that in the (exact) eq. (7-14) the 0i part is exactly the Newman- Penrose 
translation of the GHP equation which describes the 0i electromagnetic perturbation in a 
vacuum type D metric present in the paper of Lun and Fernandez©: 

[(P - P ~ P~W + 2p)-(d-r + n)(d' + 27r)]0! = , (7-15) 

which, after a tetrad transformation which sets e = 0, becomes the well known Fackerell and 
Ipser equation©. 

It is important to get the general vacuum equation for <pi in the GHP formalism, which 
can be quite easily obtanined by following the Stewart and Walker prescriptions. We only 
give the final result: 

[(p - P ~ PW - y) -(d-r- f'){& - 2r')]0! 

+ [k'p + p/t' - a'd - da' + o't + a'f' - k! 'p - /t'p]0 o (7-16) 
+ [ K p' + p' K ~~ a $ ~~ & + ra + r'a — p'k — p'k](J)2 = . 
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Translating Q7-16Q into standard Newman-Penrose formalism one gets (7-14). 

At this point the (first order) perturbation analysis can be performed. In the type D 
case, recalling that for vacuum spacetimes one has to use the test field approximation (no 
changes in the vacuum background means no changes in the tetrad vectors, spin coefficients 
and Weyl scalars), it is straightforward to show the coincidence of the linearized form of 
(7-6), (7-10) and (7-14) with the Teukolsky, Fackerell-Ipser and Lun-Fernandez equations. In 
a more general vacuum geometry the decoupling is not present and one obviously recovers 
the Stewart and Walker theorems. 



§8. Concluding remarks 

In this article we have presented a new form of the Teukolsky Master Equation and a 
new point of view for perturbations in General Relativity, starting from exact wave equa- 
tions in curved vacuum spacetimes valid for the Riemann and Maxwell tensors themselves 
and involving both the generalized and the ordinary de Rham Laplacian respectively. Once 
perturbed these equations immediately give the standard perturbation theory of Teukolsky, 
Stewart and Walker, Lousto and Campanelli, explaining also the true origin and the ge- 
ometrical meaning of the same Teukolsky Master Equation. We have discussed our exact 
wave equations both in the Newman-Penrose and Geroch-Held-Penrose formalisms, collect- 
ing most of results scattered in literature in a series of papers and books and over a period 
of more than thirty years. We have also presented a systematic study of perturbation theory 
at every order, which recovers and extends existing results on this subject. Furthermore, our 
approach can be directly generalized to nonvacuum spacetimes: 1) The de Rham Laplacian 
of the Riemann tensor then involves (in the exact regime) the Ricci tensor and curvature 
scalar terms, in turn to be re-expressed in terms of energy-momentum tensor of the sources 
or directly coupled to the de Rham equations for the other corresponding physical fields: 
Maxwell, Yang Mills, Dirac, Rarita-Schwinger0; 2) the perturbation regime can be devel- 
oped straightforwardly, recovering, for instance, once our general results are specialized to 
a vacuum type D background, the "matter perturbed sources" appearing in the works of 
Teukolsky and of Lousto and Campanelli. 

Regarding the tetrad and gauge invariance of the theory, the procedure that one must 
follow is the standard one, well discussed in the literature by Bruni, Matarrese, Mollerach 
and Sonego©'© and by© and we will adapt it to this formalism in a forthcoming paper. 

We conclude pointing out that our approach could have useful consequences for the still 
open problem of the electromagnetic and gravitational coupled perturbations of the Kerr- 
Newman spacetime^'^'^'i^ as well as for the more general problem of propagation of 
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gravitational waves on any given background. 
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Abstract 

A new version of the Teukolksy Master Equation, describing any massless field of 
different spin s = 1/2,1,3/2,2 in the Kerr black hole, is presented here in the form 
of a wave equation containing additional curvature terms. These results suggest a re- 
lation between curvature perturbation theory in general relativity and the exact wave 
equations satisfied by the Weyl and the Maxwell tensors, known in the literature as the 
de Rham-Lichnerowicz Laplacian equations. We discuss these Laplacians both in the 
Newman-Penrose formalism and in the Geroch-Held-Penrose variant for an arbitrary 
vacuum spacetime. Perturbative expansion of these wave equations results in a recur- 
sive scheme valid for higher orders. This approach, apart from the obvious implications 
for the gravitational and electromagnetic wave propagation on a curved spacetime, ex- 
plains and extends the results in the literature for perturbative analysis by clarifying 
their true origins in the exact theory. 
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1. Introduction 



In recent years there has been considerable interest in perturbations of the gravitational 
field and in other spin s massless test fields (s = 1/2, 1,3/2,2) imposed on a gravitational 
background (especially on a black hole background, which is of particular interest for study- 
ing features of gravitational collapse and consequent gravitational wave emission). The 
approaches to perturbations existing in the literature can be divided essentially into two 
classes: 1) metric perturbations and 2) curvature perturbations. 

Pioneering work in the first class was done by Regge, Wheeler and Zerilli, who produced 
important results for spherically symmetric black holes by decomposing the metric 
perturbations into multipoles which satisfy a set of (linearized) equations 4 - 1 " 9 ). However, a) 
perturbing the metric components written with respect to some coordinate system is not 
a "gauge invariant" procedure as pointed out by Stewart and Walker 10 ) and b) apart from 
spherically symmetric background solutions, this approach does not seem to be very useful. 

Working in the second class, i.e., perturbing directly the components of the curvature 
tensor, one can use frames adapted to the null principal directions of the Weyl tensor of the 
background spacetime: this allows special simplifications (even partial decoupling) of the 
equations themselves. However, this approach has problems too: perturbing the components 
of the curvature tensor with respect to some special adapted frame is not a "tetrad invariant" 
procedure and one must add this difficulty to that of the gauge invariance which is still 
present. 

The "state of the art" regarding curvature perturbations is still represented by the work 
of Teukolsky 11 ) (in the context of the Newman- Penrose formalism 12 ) , NP in the following), 
which partially received its mathematical foundation by Stewart and Walker 10 ) and was 
subjected to some extensions by other authors 13 ^ " 16 ^' 30 ) " 33 ^ . Very recently, a gauge-invariant 
second-order perturbation study for a type D vacuum geometry was initiated by Lousto and 
Campanelli 35 \ also stimulated by relevant results concerning the second order perturbations 
of a Schwarzschild black hole 36 ^ . 

We introduce a new form of the Teukolksy Master Equation, closer to a "generalized 
Klein Gordon" equation. This allows us to investigate exact wave equations on a general 
vacuum spacetime: the de Rham-Lichnerowicz equation satisfied by the Riemann tensor 
as well as the ordinary de Rham equation satisfied by the Maxwell tensor in a test field 
approximation. All the results are formulated for tensor components with respect to any 
frame and then specialized to the Newman-Penrose formalism and to its further extension, 
the GHP formalism, introduced by Geroch, Held and Penrose 33 ). Starting from the exact 
equations for the Riemann and Maxwell tensor components, it is easy to give a unified and 
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complete picture of the existing theory of perturbations to any order, framing all the known 
results (which are scattered in many different formalisms) in a more general form. In type 
D spacetimes all the known results, including those of Teukolsky , Stewart and Walker 10 - ) 
and Lousto and Campanelli 35 - ) are easily recovered and discussed. Finally, some remarks are 
given about how this approach can be generalized to half-integer spin test fields as well as 
to the case of nonvacuum spacetimes. 



52. The Newman-Penrose formalism: a brief review 



In this article we introduce the Newman-Penrose form of the generalized de Rham wave 
operator for the Weyl, Riemann and Maxwell tensors. To accomplish this, it is first nec- 
essary to specify all the formalism conventions: we will follow exactly the notations of 
Chandrasekhar 37 K For the sake of completeness we recall here some details. 

A Newman-Penrose frame is defined by four complex null vector fields ei = 1, e 2 = n, 
e 3 = m, e 4 = m*, satisfying the relations 



1 • m = 1 • m* n m n m* = 0, 
l l = n n = m m = m* ■ m* = 0, 
1 ■ n = 1, m ■ m* = —1, 

so that the components of the metric tensor in this frame are 



Vab = V 



(2-1) 
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o \ 
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-1 


Vo 





-1 


o J 



(2-2) 



With an abuse of notation (not distinguishing vectors and 1-forms) the dual frame is 

e 1 = n = e 2 , e 2 = 1 = e x , 



3 * 4 

e = — m = — e 4 , e = — m = — e 3 . 



(2-3) 



The basis vectors, thought of as directional derivative operators when acting on scalars, are 
denoted by 



e x = e 2 = D, e 2 = e 1 = A , 
e 3 = — e 4 = 5, e 4 = — e 3 = 5*, 



(2-4) 
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and the 12 complex spin coefficents are introduced through the following linear combinations 
of the 24 Ricci rotation (connection) coefficents 7 ca b = e c • V e( ,e a = —^ ac b 



K 


= 7311, 


P 


= 7314, 


e = 


2l72ii + 734i), 


a 


= 7313, 


P 


= 7243, 


7 = 




A 


= 7244, 


T 


= 7312, 


a = 


7^(7214 + 7344), 


V 


= 7242, 


7Y 


= 7241, 


P = 


2(7213 + 7343) 



As a general rule, complex conjugation of null frame components is equivalent to interchang- 
ing the indices 3 and 4. Using this prescription the inverse relations of (2-5) are 



7121 


= -(c + 0, 


7122 


= -(7 + 7*), 


7123 


= -(P + a*), 


7124 


= -(a + F), 


7131 


= -K, 


7132 


= -T, 


7133 


= -v, 


7134 


= ~P, 


7141 


= -« , 


7142 


= -r*, 


7143 


= ~P\ 


7144 


= -o- , 


7231 


= 7T , 


7232 




7233 


= X\ 


7234 


= P , 


7241 


= 7T, 


7242 




7243 


= P, 


7244 


= A, 


7341 


= (^-0, 


7342 


* 

= 7-7 , 


7343 


= 13 -a*, 


7344 


= a-(3*. 



(2-6) 



The commutation rules (Lie brackets) [e a , e b ] = C c a be c reduce to 

[A D] =(7 + 7 *)D + (e + e*)A - (r* + tt)5 - (r + 7r*)<T, 

[5, £>] =(a* +0- ti*)D + kZ\ - (p* + e - e*)5 - <7<T, 

[5, Z\] = - + (r - a* - /3)A + (// - 7 + 7*)^ + A*5*, ^ ?) 

[5*, 5] = (,/ _ p)D + (p* - p)Z\ + (a - /T)<5 + (/3 - a*)<5*. 

The 10 independent components C abcd of the Weyl tensor are represented by 5 complex 
scalars in the Newman-Penrose formalism 

% = - C 13 i 3 = - C abcd l a m b l c m d , 
*i = - C1213 = - C abcd l a n b l c m d , 

&2 = ~ C1342 = - C abcd l a m b m* c n d , (2-8) 
% = - Ci 242 = - C abcd l a n b m* c n d , 
&4 = ~ C 242 4 = - C abcd n a m* b n c m* d , 
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with the additional properties 



C1334 


— C1231 — 01, 


Cl241 


— C1443 — ?f r *, 


C1212 


=^3434 = + 


Cl234 


= (^2-^*), 


C2443 


= — C*1242 — 03, 


Cl232 


= ^2343 — — 



(2-9) 



and C1314 = C2324 = C1332 = C1442 = 0. Analogously, the 10 independent components R a b of 
the Ricci tensor are represented by the following (four real and three independent complex) 
scalars, packaged in their tracefree {&) and pure trace (A) parts 



^00 = 


- \Rii, 


022 


- ~ lR22, 


^02 = 


1 

- 2^33, 


020 


= - 2^44, 


011 = 


— ^(-^IS + R34), 


001 


= " ^13, 


0io = 


- \Ru, 


012 


= " 2^3, 


021 = 


- -R24, 


A 


24 



(2-10) 



— (^12 - R34) 



The components R a bcd of the Riemann tensor are then related to the Weyl and Ricci tensors 
by: 



-R1212 


= _ ^ 2 _ - 2<P n - 10/1, 


-Rl324 


=0 2 + 2A, 


-R1234 


=02 - 2 *, 


-^3434 


= - & 2 - 2 * + 20ii - 10/1, 


-R1313 


= -0o, 


-^2323 


^4 , 


-R1314 


= - 0oo, 


-^2324 


= - 022, 


-R3132 


=002, 


-Rl213 


= -!^i-3>oi, 


-R1334 


= -#*-$>„!, 


-Rl223 


=0 3 * + 0i2, 


-R2334 


= - 3 * - 012 , 







(2-11) 



with the additional complex conjugate relations obtained by interchanging the indices 3 and 
4. 

The six independent components F ab of the (antisymmetric) Maxwell tensor are repre- 
sented by the three complex scalars 

00 = 013 , 01 = ^(012 + 043) , 02 = 042 , (242) 
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with the inverse relations 



ab 



/ 

-<f>l - 0i 

-00 

v -<t>l 



^ + 01 



02 
02 



00 
-02 



-01 + 01 



0S \ 

-02 
01-01 

/ 



(2-13) 



Finally the Ricci and the Bianchi identities are given explicitly in the book of Chan- 
drasekhar 37 ) (section 1.8). Here they will be referred to according to the components of the 
Riemann tensor which give rise to the corresponding equation: [R a bcd] stands for the Ricci 
identity associated with R a bcd and R a b[cd\e] = for the Bianchi identities. 

Exchange of 1 <-> n, m <-> m* in the various Newman-Penrose quantities implies 



< — > ■nr dn 



D = Fd» to 



o = mro^ < — > mro^ < — > 



m*^ = 5*, 



(2-14) 



l^n m^m* 

« = 7311 < ► 7322 < ► 7422 = ~7242 = ~V, 

l^n m<->m* 

T = 7312 < ► 7321 < ► 7421 = ~7241 = TTj 

l^n m<->m* ^ 

C = 7313 < ► 7323 < > 7424 = ~7244 = "A, 

l^n m<->m* 

P = 7314 < ► 7324 < ► 7423 = "7243 = 



(2-15) 



e = ^(7211 + 734i) to i(7 122 + 7342 ) 



< > 



-(7122 + 7432) = -7, 



Oi = -(7214 + 7344) to>i(7 124 + 7344) ^(7123 + 7433) = ~P, 



(2-16) 



% = 

#2 = 



1313 < > 



UiTO? < > 



-c 



1213 



"tvl342 < > 



^2323 
^2123 
^2341 



< > 



02^41 < > 



-C2424 = $4, 

"C2124 = — Cl242 = $3) 

-^2431 = — C*1342 = $2 



(2-17) 



21 



0i = \{F l2 + F 43 ) to 1(F 21 + F 43 ) 



" ^(F 21 + F 3 4) = -0 1 . 



(2-18) 



These operations are useful since they may be used to extend a given component calcu- 
lation to other components. 
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§3. NP gravitational perturbations 

The study of perturbations of the various fields in the NP formalism is achieved by 
splitting all the relevant quantities in the form I = l A + l B + $4 = \P A + + a — 
a A + a B + ...,D = D A + D B + etc., where the A terms are the background and the B's are 
small perturbations etc. The full set of perturbative equations is obtained inserting first of 
all these splitted quanties in the basic equations of the theory (Ricci and Bianchi identities, 
Maxwell, Dirac, Rarita-Schwinger equations etc..) and keeping only first order terms. After 
some "ad hoc" algebraic manipulations, which generate second order differential relations, 
one usually obtains second order coupled linear PDE for the curvature quantities. However 
for some very special spacetimes, as for Petrov type D ones (i.e. Kerr metric), some PDE's 
are decoulpled. In the following we'll present an example of this procedure, applied to a 
generic vacuum type D metric, for which it follows that: 

$ A = & A = = & A = (3-1) 

K A = a A = v A = X A = . (3-2) 

We start from the exact NP relations: u ); 

(5* - Aa + ir)% -(D-Ap- 2e)#i - 3^ 2 = (5 + n* - 2a* - 2/3)<P 00 

—(D - 2e - 2p*)<P 01 + 2a<P 10 - 2k<P u - k*<P 02 (3-3) 

(A - 4 7 + p)% -(5 -At- 2(3)$! - 3a& 2 = (S + 2tt* - 2{3)<P i 

-{D -2e + 2e* - p*)<P 02 - \*<P 00 + 2a$ u - 2k$ 12 (3-4) 



(D - p - p* - 3e + e*V - (5 - r + n* - a* - 3(3)k -% = . (3-5) 

In the tetrad language, (3-3) represents i?i3[i3|4] = 0, (3-4) represents i?i3[i3|2] = (they 
are two components of the Bianchi identities) and (3-5) is the Ricci identity component 
-R1313 = 0. In this notation we have: 

OO = —R^ Z<T = AixT^r ee AnTu (3-6) 

and in the same sense one has to look at the other fyj. Obviously, in (3-6) n is not a spin 
coefficient but it is the usual mathematical constant coming from the right hand side of 
the Einstein equations. Because we are in vacuo, after the perturbative expansion, in these 
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relations we have that \P A , \P A , a A , n A , A and all the <& A m are zero and consequently we 
get: 

(5* - Aa + n) A V B -(D-Ap- 2e) A V B - 3k b V a = 
An[(6 + vr* - 2a* - 2(3) A T B - (D - 2e - 2p*) A T* ] 



(3-7) 



(A - 4 7 + - (<J - 4r - 2/5) A ^f - 3a B ^ = 

4tt[(5 + 2tt* - 2/3) - (D - 2e + 2e* - p*) A T r B ' 



(3-8) 



mmJ 



(D - p - p* - 3e + e*) A cr B - (5 - r + n* - a* - 3(3) a k b - V B = . (3-9) 

In the following, we'll obmit the A superscript for the background quantities to simplify the 
notation. The type D background metric satisfy the relations: 

D^ 2 = 3p^ 2 (3-10) 

6& 2 = 3t& 2 , (3-11) 
so multiplying (3-9) by \P 2 and taking into account that 

V 2 Da B = D(V 2 a B ) - a B D& 2 = D(V 2 a B ) - a B 3p& 2 (3-12) 
& 2 5 K B = 5(V 2 k b ) - K B m 2 = 5(V 2 k b ) - k b 3t& 2 (3-13) 
we obtain the relation: 

(D — 3e + e* — Ap — p*)V 2 a B - (5 + tt* - a* - 3(3 - At)V 2 k b - V B V 2 = 0. (3-14) 

We have now to eliminate \P B from the eqns (3-7) and (3-8). To accomplish this task we use 
the commuting relation found by Teukolsky and valid for any type D metric: 

[D-(p+ l)e + e* + qp - p*] (5 - p(3 + qr) 

-[6 - (p + l)(3 - a* + tt* + qr](D - pe + qp) = (3-15) 

where p and q are two generic constants. We apply now (D — 3e + e* — 4p — p*) on (3-8) and 
(5 + it* — a* — 3(3 — At) on (3-7) and we subtract the two resulting equations. We can now 
use (3-15) with the choice p = 2 and q = —4 to eliminate the \P± term. The combination o B 
and n B remained is exactly the same present in eq. (3-14), and consequently these quantities 
can be eliminate in favor of the $2$cf term. The final equation is: 

[(D - 3e + e* - Ap - p*)(A + p - 4 7 ) - (5 + vr* - a* - 3(3 - At) 

x (5* + tt - Aa) - W 2 ]$o = 4ttT (3-16) 
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where we have posed: 

T = {5 + 7i* - a* - 3(3 - At) [(£> - 2e - 2p*)T* - (5 + tt* - 2a* - 2/3)7}f ] 
+ (£, _ 3e + e * _ 4p - p *)[(<5 + 2tt* - 20)7* -(D-2e + 2e* - p*)T* m ] . 

(3-17) 

In the same way one can build the other relations to perform the perturbative analisys 
of a type D vacuum spacetime. We point out that the same technique can be applied to 
the NP Maxwell equations to decouple them as well to the neutrino and Rarita-Schwinger 
ones 11 )' 25 ) It is important to stress that this perturbative theory has a purely "algebraic and 
manipulative" nature, without taking into account the geometry underlying the problem. 
In the following we'll give to this subject a new geometrical meaning, framing it into the 
context of the field theory. But to accomplish this task we have to adapt the previous studies 
to the Kerr solution, which describes a rotating black hole of specific angular momentum a 
and mass M, with the condition of flat spacetime at infinity. 

§4. The Teukolsky Master Equation 

The Kerr solution in Boyer-Lindquist coordinates is given by 21 -*: 

, , / 2Mr\ , 2 4aMrsin 2 6> , , , E , „ 
ds 2 = f 1 J^) dt + ^ dtd( ^ ~ -£ dr 

a 2 sin 2 9 



- Ed9 2 - 2Mr 



r 2 + a 2 + 



sin 2 9d(f) 2 (4-1) 



where as usual: 



A = r 2 - 2Mr + a 2 , E = r 2 + a 2 cos 2 9 (4-2) 



For any Petrov type D metric, and in particular for the Kerr solution, simplifications occur 
working in a tetrad adapted to the two repeated principal null directions of the corresponding 
Weyl tensor. In the NP formalism the following quantities 1? ) completely describe the Kerr 
solution (in this section we use the label A over some quantities in view of a perturbative 
analysis, as will be clear in the following). The Kinnersley tetrad 18 ^: 

(n A = ^[r 2 + a 2 ,A,0,a] 

(nn A =^[r 2 + a 2 ,-A,0,a] (4-3) 
{m") A =^ [ia sin 0,0,1, ' 



y/2(r + ia cos 9) ' sin 9 
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with the 4th tetrad vector being (m* M ) A , the complex conjugate of [m^) A . The Weyl tensor 
is represented by: 

g/ A = = = & A = 

(44) 

= M{p A f 

while the Ricci tensor and the curvature scalar are identically zero because of the vacuum 
spacetime. The spin coefficents are given by: 



k = a — a = v — e = , 



A — 1 A —iap A p* A sin 9 

^ (r — mcosfl)' ^[2 

A -p* A cot9 A ia(p A ) 2 sin9 

P = n = V2 ' (4 ' 5) 

, (p A ) 2 p* A A A _ A p A p* A (r - M) 

V - 2 7 -/i + 2 



a A = tx a - (3* A 



It is a relevant result that all the massless (non trivial) perturbations of different spin of 
the Kerr metric are included into the Teukolsky Master Equation (TME) n ) which has the 
following form (ir is the mathematical constant): 



2\2 



(r 2 + a 2 ) 



-A~ s 



A 
d_ 

Or 



a 2 sin 2 9 



d 2 ib AMar d 2 ib 

+ a + 



\s+l 



dip" 
dr 



-2s 



~M(r 2 - a 2 ) 



dt 2 
1 d 
sin 9 89 



r — ia cos 9 



A 

sin 9 



dtdcp 
dip" 
~d9 



a 
A 



2s 



sin 2 9 
a(r - M) 



d 2 j) 
i cos 9 

+ 



A sin 2 9 

^ + (s 2 cot 2 fl- S W = 47rrT. 
at v ' 



dip 
80 



(4-6) 



The following table shows the various quantities involved in the Teukolsky Master Equa- 
tion (T a and J a are respectively gravitational and electromagnetic sources): 



TME quantities 
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i> 


S 


T 




<p 





see ref. 20 ^ 




Xo 


+ 1/2 


see ref. 34 ^ 


p 




-1/2 


see ref. 34 ) 




00 


+ 1 


Jo 


p 




- 1 


P~ 2 -h 




0> 


+ 3/2 


see ref. 25 ^ 


p~ 




-3/2 


see ref. 25 ** 




% 


+ 2 


2T 


p 


~ 4 % 


- 2 


2p~ 4 T 4 



(4-7) 



Table 1: The physical field component ip, spin weight s and source terms T for eq. (4-6) 

It is possible to prove that the other NP quantities not described by this equation are 
related to trivial perturbations in the Kerr background 30 )' 4 % and consequently the study 
of the Teukolsky Master Equation alone allows practically the recovery of all the physical 
quantities. The reader can easily verify that eq. (3T6) coincides with (4-6) for s = 2. It is 
remarkable the separability of this equation in the form: 



^(f , r, 8, 0) = e- iuJt e im,i, R(r)Y(9), 



(4- 



by using the so called "spin weighted spheroidal harmonics" which, as it is well known from 
molecular physics 19 ^ , satisfy the following spectral problem: 



(Ho + Hi)Y(9) = -EY{9) 



(4-9) 



with 



Hn 



1 d ( . d_ 



m 2 + s 2 + 2ms cos 9 
shi 2 ! 



(4-10) 



Hi = a 2 uj 2 cos 1 9 — 2aujs cos 9, 



(4-11) 



and E = E 8t i >m and Y = Y s> i >m are the angular eigenvalues and eigenfunctions which in 
general can only be obtained perturbatively. At this point, the perturbation of the Kerr 
black hole has become a typical Quantum Mechanics work in which one has to solve a 
Sturm-Liouville problem for two coupled radial and angular ODE's. The Teukolsky Master 
Equation can be cast in a more compact form by introducing a "connection vector" whose 
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components are: 



1 

E 



M(r 2 - a 2 ) 



(r + ia cos 6) 



jpr 

r e = o 



-(r-M) 



= -— 



i 

E 



a(r — M) . cos 9 
A + % ^9 



It's easy to prove that: 



vr u = -- 



i 



i 



r>*r u = ^ cot 2 e + 4^ 



and consequently the Teukolsky Master Equation assumes the form: 



[(V M + si^)(V M + sr M ) - 4sV 2 A ]^ (s) = 4vrT 



(4-12) 



(4-13) 



(4-14) 



where is the only non vanishing NP component of the Weyl tensor in the Kerr background 
in the Kinnersley tetrad (4-4). Equation (4-14) gives a common structure for these massless 
fields in the Kerr background variating the "s" index. In fact, the first part in the lhs 
represents (formally) a D'Alembertian, corrected by taking into account the spin-weight and 
the second one is a curvature (Weyl) term linked to the "s" index too. This particular form 
of the Teukolsky Master Equation forces us to extend this analysis in the next sections. 



§5. Generalized wave equations 

The compact form of Teukolsky Master Equation we have shown suggests a connec- 
tion between the perturbation theory and a sort of generalized wave equations which differ 
from the standard ones by curvature terms. In fact generalized wave operators are know 
as De Rham-Lichnerowicz Laplacians and the curvature terms which make them different 
from the ordinary ones are given by the Weitzenbock formulas. Mostly known examples in 
electromagnetism are 

• the wave equation for the vector potential 21 ) : 

V a VVL M - R^A X = -AnJ^ , V a A a = (5-1) 

• the wave equation for the Maxwell tensor 24 ^: 

V%F„ A + R^xF"* - R\F vp + R%F Xp = -8ttV [m J v] (5-2) 
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while for the gravitational case one has 



• the wave equation for the metric perturbations (gravitational waves) 



21). 



2R a ( t j,h l/ ) a = 0, 



• the wave equation for the Riemann Tensor 21 ) (in vacuum for simplicity) 




(5-3) 



V M V M -R a/ 3 7 5 — RappaR^s 



pa 




) = 0. 



(5-4) 



These equations are "non minimal," in the sense that they cannot be recovered by a minimal 
substitution from their flat space counterparts. A similar situation holds in the standard 
Quantum Field Theory for the electromagnetic Dirac equation. In fact, applying for instance 
to the Dirac equation an "ad hoc" first order differential operator one gets the second order 
Dirac equation 27 ^ 28 ) 



where the notation is obvious. It is easy to recognize in eq. (5-5) a generalized Laplacian and 
a curvature (Maxwell) term applied to the spinor. Moreover this equation is "non minimal" , 
in the sense that the curvature (Maxwell) term cannot be recoverered by electromagnetic 
minimal substitution in the standard Klein-Gordon equation for the spinor components. 
The analogous second order Dirac equation in presence of a gravitational field also has a non 
minimal curvature term 22 )' 23 '' and reduces to the form: 



Quite recently it was shown by Mohanty and Prasanna 24 ), that the Maxwell tensor wave 
equation (5-2) could be framed into the context of the QED vacuum polarization in a back- 
ground gravitational field. They start from a classical work of Drummond and Hathrell 26 ^ 
which derives higher derivative coupling arising from the QED loop corrections to the 
graviton-photon vertex, encoded in the Lagrangian, valid for u 2 < a/m 2 e (a is the fine 
structure constant and m e is the electron mass): 



with the coefficients a = — , b = ffjjf, c = — yfgf • The physics of these "tidal" terms 
is that the effect of virtual electron loops gives to the photon a typical size proportional to 





(5-6) 




-F^ + aRF^ + bR^F^F\ + cR^F^F 



(5-7) 
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the Compton wavelenght of the electron A c . The field equations obtained from the modified 
lagrangian (5-7) are: 

V M F^ + W = (5-8) 

where Q v represents the non linear quantum correction the explicit (long) form of which is 
given in 24 ^. Manipulating the previous relation one gets the wave equation: 

VV,^ + {R P ^FW + R\F vp - R p v F Xp \ + [V V Q X - V X Q V ] = . (5-9) 

In the range of validity of the effective action the Riemann and Ricci terms in the curl 
brackets of (5-9) which are present in Einstein's gravity are larger than the terms in the 
square brackets, arising from the loop corrections, and consequently it is clear that the 
classical second order wave equation for the Maxwell tensor describes with good accuracy 
this quantum phenomenon too. 

§6. Vacuum gravitational fields 
The Weyl and the Riemann tensor are connected by the relation: 

Ral3~f8 = Cal3yS + -(ga^R/3S ~ 9^yR a S ~ 9asR/3j + QpsRa-f) 

-^{9ay9l3S ~ 9aS9lh) R , (64) 

In a vacuum spacetime, R a p = 0, R = and consequently the Riemann tensor coincides with 
the Weyl one. The Weyl tensor in vacuum satisfies the generalized de Rham wave equation 
(the vertical bar indicates the covariant derivative) 

(6-2) 

The name "generalized de Rham wave equation" comes from the nice compact form which 
can be given to this equation 

AdR) C ai a 2 a 3 a 4 = ($D + D5)C aia2a3a4 = , (6-3) 

by using the divergence operator 5 and the covariant exterior derivative D 39 ) " 41 ) (the diver- 
gence operator S and the covariant exterior derivative D must not be confused with the NP 
frame derivatives denoted by the same symbols). This is the most natural generalization to 
tensor-valued forms of the ordinary de Rham wave operator 

A (dR) = 5d + d6, (6-4) 
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whose action is denned only on ordinary differential forms. A detailed study of the general- 
ized de Rham operator will be discussed elsewhere 38 - ) . 

The explicit form taken by [Z\(dR)C]a6cd = in the Newman-Penrose formalism is now 
derived and simplified somewhat using the commutator identities to eliminate some second 
derivatives and the Bianchi and Ricci identities to eliminate some coupling terms between 
different Weyl tensor scalars. 



AdR)C]i3i3 = (wave equation for %) 



(6-5) 



Recalling that $b = — C1313, select the 1313 component of eq. (6-2) 

[—DA -AD + 5*5 + 55* + (9 7 + 7* - p - p*)D + (7e - e* + p* + p)A 
+ ((3* + n -9a- t*)5 + (tt* - r - 7(3 - a*)5* + 4(£> 7 + Ae - 5a - 5*(3) 
+ 4(aa* - p*7 + t*(3 + r7r - f3*f3 - 8^e - £7* + 7 e* - an* + pe - f3n - a\ 
+ Ku + ra + 8a (3 - pp - p 7 + p*e + 3& 2 )}% + 4[2p5 + 2a5* - 2rD - 2k A 
+ 5p — Dr + 5* a — Ak + n*p — pa* + 5re — re* — hpj3 + rp* + an — up* 
-Kfi- 7aa + /t 7 * + 7^7 - or* + a[3* - + [2A(ap - kt)}& 2 = . 

This equation only contains second order directional derivatives of $0, which can be rewritten 
first using the identities 

DA + AD = 2D A + [A, D] 

55* + 5*5 = 255* + {5*, 5} (6-6) 

and then the commutators in these expressions can be replaced by first order derivatives 
using the commutator identities (2-7)i and (2-7) 4 , in order to reduce the number of second 
order differential operators acting on %. 
Equation (6-5) then becomes 

[-2DA + 255* + (87 - 2p)D + (6e - 2e* + 2p*)Z\ + (2tt - 8a)5 
+ (2tt* - 6/3 - 2a*)5* + 4(£> 7 + Ae - 5a - 5*13) + 4(aa* - p* 7 
+ T*i3 + T7T - 13*13 - 87c - £7* + 7e* - air* + ae - 3n - aX 

(6-7) 

+ ku + ra + 8aj3 - pp - p>y + p*e + 3^ 2 )]^o + 4[2p5 + 2a5* - 2rD - 2kA 
+ 5p — Dr + 5* a — Ak + n* p — pa* + 5re — re* — 5p/3 + rp* + an — Kp* 
- Kp- 7aa + K>y* + 7k^ - ar* + af3* - 3<?i]!Ji + [24(crp - kt)\P 2 = 

Next the Bianchi identities can be used to simplify this expression and introduce terms 
involving % to decouple it whenever possible. This task is achieved by the elimination of 
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two directional derivatives of #1 in eq. (6-7) using two Bianchi identities 

#i3[i3|4] = : D&! - 5*% + (4a - tt)% - 2(2p + e)^i + 3k& 2 = , 
i2i3[i3|2] = : 6$! - A% + (4 7 - n)V - 2(2r + (3)^ + 3a& 2 = . (6-8) 

Equation (6-7) then becomes 

[-2£>Z\ + 255* + 2(4 7 - p)D + 2(3e - e* + p* + 4p)Z\ + 2(tt - 4a)<5 
+ 2(tt* - 3(3 - a* - 4r)<T + 4(£> 7 + Z\e - 5a - 5*13) + 4(aa* - p* 7 
+ r*/3 - T7T - tf*/? - 8 7 e - e 7 * + 7 e* - avr* + pe - (3tx - aX 

(6-9) 

+ kv + 9ra + 8a/9 + pp - 9p 7 + p*e + W 2 )]% + 4[2cr5* - 2/tZ\ 

+ 5p — Dt + 5*o" — Ak + 7r*p — pa* + re — re* — p/3 + rp* + o"7r — /tp* 

- Kfx- 7cra + /t 7 * + 7/t 7 - or* + a"/?* - 3«Pi]«?i = . 

Finally this (exact) equation can be given a Teukolsky-like form which is very helpful in 
type D geometries where most of the known results have been derived. To accomplish this 
additional step, multiply the equation by —1/2 and use three Ricci identities in two different 
ways: a) in (6-9) first use the identities 

[l/2(#i2i 2 - #3412)] : At = D 1 - a(r + tt*) - (3{r* + vr) + 

+ 7 (e + e*) + e( 7 + 7*) - rvr + uk - & 2 , 

6-10 

[1/2(^1234 - #3434)] : S*(3 = 5a- (pp - \a) - aa* - (3(3*+ 

+ 2a(3 - 7 (p - p*) - e(p - p*) + V 2 ; 
b) rewrite the other Ricci identity [#2431] in the following form 
[#2431] : Dfi-Sir- (p*p + o-A) - 71(71* - a* + (3) + p(e + e*) + uk - & 2 = , (6-11) 

multiply it by \P and c) add this (which is actually zero) to the wave equation. The result 
is 

[(D - 3e + e* - 4p - p*)(A + p - 4 7 ) - (5 + tt* - a* - 3(3 - 4r) 
(5* + n- 4a) - 3W 2 }% + 3(<t A - kv)% + 2[2kA - 2o5* - 5*a + An 

(6-12) 

— 5p + Dt — 7/c 7 — (3* a + 7aa + r*cr + up* + up — pix* + re* + (3p 

- er - k 7 * - Tin - rp* + pa* + Wi\^i = ; 
d) finally use in eq. (6-12) two other Ricci identities 

[#3143] : Sp - 5* a = p(a* + (3)- a(3a - (3*) + r(p - p*) + n(p - p*) - ^ , 
[#1312] : Dt-Ak = P (t + tt*) + a(r* + tt) + r(e - e*) - «(3 7 + 7*) + 3^. (6-13) 
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The final form for eq. (6-12) is then 

[(D - 3e + e* - 4p - p*)(A + fj, - 4 7 ) - (5 + vr* - a* - 3/3 - At) 

(5* + n- 4a) - 3V 2 ]&o + 3(aA - ku)% + [4kA - 4a5* - 45*a + 4Ak (6-14) 

- 20/t 7 - 4/3 V + 20<ra + 4tV + 4k//* - 4/t 7 * + lO^i = , 

which is the exact form for the same equation for #o found by Teukolsky in the very special 
case of a type D geometry in the first order perturbation regime (see below). Moreover 
this equation coincides with that obtained by Stewart and Walker w > following a different 
approach and in the Geroch-Held- Penrose (GHP) formalism 33 ) 

( p'p - &'d - p p - 5pp' + f d + hrd' + 4aa' - 4kk' - 10V 2 )% 

+ [4£'/c - 4d'a - 4(p' - 2p')« + 4(f - 2r> + lO^i (6-15) 

+ (-4(7p + 4k<9 - 12/tr + 12pa)^ 2 = . 

The equivalence of (6-14) and (6-15) is now explicitly shown. Recall the action of the 
two-index {p, g}-operators p, p', d, & on two-index {p, g}-objects. Taking into account the 
weights of the !^ n 's (i.e., {p, q} = {4 — 2n, 0}), those of the spin coefficients and the definition 
of the weighted derivative operator of an object rj (which increases its {p, q} indices) 



pr? = (D — pe — qe*)r), increasing of type {+1, +1}, 

\>'r] = (A — p'j — q'j*)'!], increasing of type { — 1, —1}, 

drj = (5 — p(3 — qa*)r), increasing of type {+1, —1}, 

d'rj = (5* — pa — q(3*)r), increasing of type { — 1, +1}, 

eq. (6-15) becomes: 

[(A - 57 - 7*)(£> - 4e) - (5* - ha + (3*)(5 - 4/3) + p*(D - 4e) 

- 5p(A - 4 7 ) +r*(5- 4/3) + 5r(5* - 4a) - 4(<rA - kv) - 1W 2 }% 
+ (-4<tD + 4k5 - YIkt + \2pa)^ 2 + [4(Z\ - 3 7 - 7 *)k 

- 4(5* -3a + (3*)a - 4(-/i* + 2/i)/t + 4(r* + 2n)a + 10^)]^ = , 



(6-16) 



(6-17) 



where the second order operator has to be rewritten via the commutation rules (2-7). 

Then, by using the two Ricci identities 1/2[R 1212 — -R3412] and 1/2[R 12U — -R3434] (see 
eq. (6-10)), adding equation [-R2431] multiplied by \P and using the two Bianchi identities 
#i3[43|2] = and i?i 3 [2i|4] = 0, namely: 

A&! - 6& 2 - IS&o ~ 2(7 - + 3r^ 2 - 2<J& 3 = , #13[43|2]=0, 

(6-18) 

5*^i - DV 2 - X% + 2(vr - a)&! + 3p^ 2 - 2/t^ 3 = , #i 3 [2i|4] = , 
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one gets exactly eq. (6-14). 



[AdR)C]2424 = 



(wave equation for ^4) 
The equation for $4 can be immediately obtained from the equation for ^ using the ex- 
change symmetry 1 <-> n, m <-> m* of the Newman- Penrose quantities. Under this operation, 
eq. (6-14) becomes: 

[(A + 3 7 - 7* + 4/i + p*)(D + 4e - p) - (5* - r* + 0* + 3a + An) 

(5 - r + 4/?) - 3V 2 ]V A + 3(aA - ki/)^ 4 + [4A<5 - 4zaD + A5X - ADv (6-19) 

- 20z/e - 4a*A + 20A/3 + 4tt*A + Avp* - Ave* + 10# 3 ]#3 = . 



Working in the GHP formalism, the corresponding relation found by Stewart and Walker 



10 ) is: 



[( p'p - d'd - (Ap' + p')p - pV + (4r' + f)d + r& + App' - Att' - 2^ 2 ]^ 4 

+ [4p«/ - Ada' - 4(p - 2 P )k' + 4(f - 2r)a' + lOlPj,]?^ (6-20) 

+ [-4cr / p / + Ak'& - 12kV + 12pV]^ 2 = 0. 



[^(dR)C; 



1213 



= 



(wave equation for &1) 



Selecting the component C1213 ($1 = — C1213) in eq. (6-2) one has 



[55* + 5*5 - DA - AD + (7* + 5 7 - p - p*)D + (p + p* + 3e - e*)A 
+ (tt + (3* - t* - 5a)5 + (tt* - 3/5 - r - a*)<T + 2(L>7 + Z\e - 5a - <S*/3 
+ 5r7r + ft«* — 47c — p7 — p*7 + p*e — €7* — + 7e* — /?7r 

- 5a A + /5r* - a7r* + pe + 4a/3 - 5pp + ra + 5/tz/ - 3^ 2 )]^i 

+ [2vD + 2ttZ\ - 2X5 - 2p5* + Dv + An - 5X - 5*p + r*p - z/p* + tX (6-21) 
+ vrp* + Aa* - 777T - 5ve + 7ap - 7r7* + ve* - vp + 5Xf3 - Att* - p(3* + 6^3] $0 
+ 3[-2rD - 2kA + 2p5 + 2<7<T - Dr - Ak + 5p + 5*a + /3V - /3p - ar* 

— /tp* + K7* — 3<ra — re* + er — a*p + rp* + 37K — p/t + o"7r + p7r*]$2 
+ 12(ap - kt)V 3 = . 

Using the commutation rules (2-7), multiplying by —1/2, and using the six Bianchi 

identities i?l3[13|4] = 0, -Ri3[i3|2] = 0, i?l3[21|4] = 0, i?l3[43|2] = 0, i?42[13|2] = i?42[13|4] = 

and the five Ricci identities [-R3143], [-R1312], [-R2443], [-R2421], [K-R1234 - -R3434))], eq. (6-21) 
becomes: 
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[DA - 55* + 2(p - i)D + (-3p - p* - e + e*)A + 2(a - tt)5 + ([3 + a* + 3r - n*)5* 

— Dj — Ae + 25a + 2p*7 — ^e* + cot* + kv + e7* + 6p7 + 47c — 6pp — 7ra + 5t7t — /5r* 

- 2aa* - 2a(3 + 3/3?r - 4pe + 7# 2 ]#i + [A5 - i/D + 5\ - Dv - 3rA + ve - a*\ (6-22) 
+ 3pz/ - ve* + Att* + up* -p\- 2V 3 }% + 3[/t5 - oD + 2ap + 2re/3 - 2kt - 2ae]^ 3 

+ 3[Ak - 5*a - 2pK + 2oix - a[3* + np* - 3«7 + or* - k^* + 3a(T]^ 2 = 0. 

Now adding to the previous relation the equation [-R2431] (see eq. (6-11)) multiplied 2$i and 
using another Ricci identity, [§(#1212 — -R3412))], (6-22) takes the compact form: 



[(£> — e + e* — 3p — p*){A + 2p - 2 7 ) - {5 - a* - + n* - 3r)(5* + 2vr - 2a) 
+ 6^ 2 ]^i + 2(ki/ - (tA)^i + [A5 - vD + 5\ - Dv - 3rA + i/e - a*A + 3pz/ 
- z/e* + Att* + z/p* - (3X - 2V 3 }% + 3[k5 - aD + 2<rp + 2k(3 - 2kt - 2ae}& 3 
+ 3[Ak - 5*a - 2pK + 2air - a(3* + Kp* - 3^7 + or* - K7* + 3aa]& 2 = . 
The analogous equation in the GHP formulation is 

[(p - p - 3p)(p' - 2p) -(d-T - 3t){& - 2t') + 2{aa' - kk') + 6V 2 ]&i 
+ [3k& - 3(jp + Qap - Qkt]V 3 

+ [pre' + «/p - <j'd - + cr'f ' + 3a t - k'p - 3k' p - 2& 3 }V 
+ [3j>'/c - 35V + 6/tp' - 6ar' + 3cxf - 3/tp']^ 2 = . 



(6-23) 



(6-24) 



AdR)C]l242 = 



(wave equation for 
$3 can be easily obtained, as was done for ^4, applying the 



The equation for C1242 = 
interchange symmetry 1 <-> n, m <-> m* to (6-23), yielding: 

[(A + 7 - 7* + 3p + //)(£> - 2p + 2e) - (5* + /T + a - r* + 3n)(5 - 2r + 2(3) 
+ 6^ 2 ]^3 + 2(kv - a\)& 3 + [kA - a5* + Ak- 5*a - 3aiT + K>y - (3*a + 3p,K 
- kj* + ar* + np* -aa- 2^i]^ 4 + 3[AZ\ - v5* + 2\p + 2va - 2vix - 2\^\P\ 
+ 3\5\ -Dv- 2pv + 2Ar - \a* + vp* - 3ve + Att* - ve* + 3(3\}^ 2 = . 
The analogous equation in the GHP formulation is 



(6-25) 



[(p' - p - 3p')(p - 2p) - (& - f - 3t'){8 - 2t) + 2(a'a - k'k) + 6& 2 }W 3 
+ [3k 1 d' - 3<7'p' + 6(7'/ - 6kV]!^i 

+ [p'/t + «!>' - 0-9' - a'cr + erf + 3(7t' - /tp' - 3/tp' - 2$i]^ 4 
+ [3f>/c' - 35a' + 6k' p - 6<tV + 3cr'f ' - 3/c'p]!^ = . 



(6-26) 
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[Z\( dR )C]i342 = (wave equation for W 2 ) 



Selecting the component 1342 of (6-2) (\P 2 = _ C1342) one has 

[—DA -AD + 55* + 5*5 + (7 + 7* - p - p*)D + (p + p* - e - e*)A 

+ (n + f3* -a- r*)5 + (tt* + (3 - a* - t)5* + YIttx - Ylpp + YLvk - 12Xa - W 2 \£ 2 

+ [AvD + An A - AX5 - Ap5* + 2Du + 2 An - 25X - 25* p + 2pr* - 6jn + 2ue* - 2p*u 

- 2pp* + 2X(3 + 2\a* + 2np* - 2n-f* - 2vt + 2rA - 2pv - 2\n* + 6ap + 4^]^ (6-27) 
+ [-ArD - AkA + Ap5 + Aa5* - 2L>r - 2Ak + 25p + 25* a + 2tt*p - 2e*r + 2rp* 

- 2/i/t + 2aa + 6p/3 - 2rV + 2a [3* + 2na - 2np* - 2«7 - 6er - 2pa* + 2/t 7 *]^ 3 
+ [Apa - Atk]& a + [AXp - Avn + 2& 4 ]& Q = . 

Using the commutation rules (2-7), multiplying (6-27) by —1/2, and using the four Bianchi 
identities -Ri 3 [2i|4] = 0, i?i 3 [43|2] = 0, -R 42 [i3|2] = -R 42 [i3|4] = 0, the four Ricci identities [-R3143], 
[-R1312], [-R2443], [-R2421] and adding eq. (6-11) once multiplied by 3$2, relation (6-27) takes 
the form: 

[(D + e + e* - 2p - p*){A + 3p) - (5 + n* - 2r + (3 - a*){5* + 3n)}V 2 
+ (3Xa - 3ukW 2 + 2\-vD + X5-Du + 5X + up* + Xn* - 2rX 

(6-28) 

- Xa* + XI3- ue* -ue + 2pv + iPj,]?^ + 2[kA - a5* + Ak 

- 5*a + aa + up* + 2pn + r*a - K'j* - - 2na - a{3*]& 3 - = 

The corresponding equation in the GHP formalism is 

[(b -2p- p)(b' - 3p') -{d-2r- f')(d' - 3r') + 3(kk' - aa')]V 2 

+2[k'\> + b/?' - a'd - da' + 2ra' + a'f' - 2 P k' - n'p}^ (6-29) 
+2[reb' + b'/t - d'a - ad' + ^1 + 2r' a - 2p'n + fa - Kp']% - = . 

At this point the complete set of (exact) wave equations satisfied by the Weyl scalars has 
been given. We can now perform their perturbation expansion. 

6.1. The perturbation theory 
Introducing the operators: 

t = [(D - 3e + e* - Ap - p*)(A + p - 4 7 ) - (5 + n* - a* - 3(3 - At)(5* + n - 4a) - 3^ 2 ], 
W = 3(aX- ku), (6-30) 
Vi = [AkA - Aa5* - A5*a + A An - 20/t7 - A/3* a + 20cra + At* a + Anp* - 4/t 7 * + 10#i], 
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where the letter % here stands for the operator on \Pq found by Teukolsky in the perturbation 
regime; this is a notation extending that due to Lousto and Campanelli 35 -* . In the same way 
a generic operator x acting on will be denoted by x%- Thus, eq. (6-14) can be written as: 

(f + W )% + = . (6-31) 

Performing the standard perturbation expansion of the various quantities a la Teukolsky: 

& Q = $W + + ^(2) + ^ a = a (0) + a (l) + ; D = D (0) + £,(1) + ... and so Qn> one finds 

[(T (0) + f (1) + f (2) + ...) + (wf + yp« + + ...)](^ (0) + ^o (1) + + ■■■) 

+ (^ 0) + + + ...)(^ (0) + + ^ (2) + ...) = o , (6 ' 32) 

giving an inductive procedure for obtaining a hierarchy of equations of the following type: 

(f (0) + wf)^ 0) + A (0 M 0) = o, 

(f (1) + wi 1 V (0) + (f (0) + wi°V (1) + (A (0) ^i (1) + v{M 0) ) = o, 



(6-33) 



(f (n) + wi n V (0) + (T (0) + H#V (n) + (Pi n) ^i (0) + Pi 0) ^ (n) ) = s£ n) , 



where 



n— 1 

= - £ [(f (n - p) + W,S n - p V (p) + pJ n -%^] , (n > 2) . (6-34) 
P =i 

For every n > 2 perturbation equation, the "source term" Sq^ is composed of quantities at 
most of order (n — 1), which in principle are known. 

For a type D background geometry it is well known that 



= <t(°) = A<°> = !/(°) = 0, 



(6-35) 



so VVq - 1 = Wq 1 - 1 = "Pi " 1 = and the recurrence scheme (6-33) becomes: 

= 0, 

T (0 ^o (1) = 0, 



(6-36) 



<f"( ),7/( n ) _ c( n ) 
1q — Oq 
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One can study the perturbation of eq. (6-19) exactly in the same way as before. In fact, 
defining the operators 



% = [(A + 3 7 - 7* + 4/i + p*){D + 4e - p) - (5* - r* + /F + 3a + 4vr) 
(8 - t + 40) - 3V 2 ], 

~ (6-37) 
V 3 = [4X5 - 4vD + 45 A - 4Dv - 20ue - 4a* A + 20A/3 + 4vr*A + 4up* - 4ue* + 10# 3 ], 

>V 4 = 3(crA - ku), 

(where the letter 71 here stands for the operator on $4 found by Teukolsky), eq. (6-19) can 
be written as 

(t 4 + W 4 )^4 + V 3 & 3 = . (6-38) 
The hierarchy of equations in this case is: 

(fi 0) + wi 0) )^ 0) + pf^ 0) = o, 

(t 4 (i) + wivf + (t 4 (o) + wi°y 4 (i) + (pf^ 3 (1) + = 0, 

(6-39) 

(f 4 (n) + vPi n Vi 0) + & 0) + vPi Vi n) + (Pfty 0) + ^ 0) ^ (n) ) = $i n \ 

where 

n— 1 

5| n) = - £ [(T 4 (n ~ p) + Wi n ~ p V 4 (p) + pf^Vf] , (n > 2) . (6-40) 
P =i 

In the case of a type D geometry w|°' ) = = = and the set of equations (6-39) 

becomes: 



= 0, 



(6-41) 



The previous results are an alternative but equivalent formulation of the work of Lousto and 
Campanelli concerning the higher order perturbations of a type D vacuum metric. The same 
perturbation formulation can easily be written for eqs. (6-23), (6-25) and (6-28) allowing a 
systematic formulation of the gravitational perturbation for every vacuum spacetime. In 
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fact by introducing the operators: 

Ai = [(D — e + e* — 3p — p*)(A + 2p - 2 7 ) - (5 - a* - (3 + n* - 3t)(5* + 2ix - 2a) + 6V 2 ], 
§! = 2(Ku-a\), 

K = [X5 -vD + 5X- Dv - 3rA + ve - a*X + 3pv - ve* + Xn* + vp* - (3X - 2V 3 ], (6-42) 
Q 3 = 3[k5 -aD + 2op + 2k(3 - 2kt - 2ae], 

M 2 = 3[Ak- 5* a - 2fiK + 2otx - a (3* + k/j,* - 3«7 + or* - k^* + 3aa], 
eq. (6-23) can be rewritten as: 

(A.! + + £ ^o + Qs^s + M 2 V 2 = , (6-43) 

and one can develop a hierarchy of perturbation equations, analogously to what we did 
above. In the case of the first order perturbation of a type D background one finds 

= = = Q(°) = ^(°) = 0, (6-44) 
so that the linearization of eq. (6-43) gives 

A ( M 1] + M ( M 0) = 0, (6-45) 

which unfortunately couples Weyl scalars and perturbed spin coefficents. Relations for \J/ 3 
can be worked out in exactly the same way. For \P 2 the following operators are introduced: 

T 2 = \(D + e + e* - 2p - p*)(A + 3p) - (5 + n* - 2r + (3 - a*){5* + 3tt)], 
Q 2 = (3Xa — 3vk), 

1 (6-46) 
M = 2[-i/D + X5 - Dv + 5A + up* + Att* - 2rA - Xa* + X(3 - ve* - ve + 2pv + 3*,], 

£ 3 = 2[kZ\ - a5* + Ak- 5*a + atr + k// + + r*a - K>y* - 7/c - 2ttct - a/3*]; 
eq. (6-28) can then be rewritten as: 

(^2 + £2)^2 + + C 3 V 3 - %V 4 = . (6-47) 
Again the perturbation theory can be directly developed. In particular, for the type D case 

g 0) = §P = Afl 0) = 4 0) = , (6-48) 
and the first order perturbations become: 

T { 2 M l) + {T? + = . (6-49) 
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Thus, even in the very special case of a type D geometry, the perturbation equations for 
&i,&2 and $3 are coupled — order by order — with other perturbation quantities. This is due 
to the presence, in the exact equations, of terms involving $2 which are responsible for the 
appearance of the perturbations of the spin coefficents and of the directional derivatives. 
To obtain decoupled equations for these three Weyl scalars in the perturbation regime, in 
principle, one has to define linear operators with derivatives of order higher than the second. 

We conclude this section by pointing out that our perturbation relations (6-36) and (6-41) 
assume a more compact form if compared with those found by Lousto and Campanelli be- 
cause of the completely different approach and the number of simplifying operations (always 
explicitly specified) that we perform at the level of the exact theory. They use 'background 
operators' acting on the exact Bianchi identities and then linearize the results order by order 
(see also their remark after eq. (6) in 35 )), following the Teukolsky prescription. The way of 
using exact 'ad hoc' operators acting on the exact Bianchi identities and then linearizing the 
results was introduced by Stewart and Walker but only for first order perturbations. 

By extending such a prescription to any order one obtains our relations. The de Rham 
Laplacian mediates this task at the level of the exact theory, also explaining the nature (i.e. 
exact 'wave equations in curved spacetime' for some of the NP Weyl and Maxwell scalars) 
of the relations found algebraically by Stewart and Walker. 

§7. Vacuum spacetimes: test electromagnetic fields 

In this section we study the electromagnetic field equation in the absence of sources, 
J M = 0. This choice has to be related to the "exact" (in the sense of the field equations) 
nature of our approach, which goes beyond the perturbation level and would be corrupted 
by the introduction of a source term without the automatic appearance of a corresponding 
term in the Riemann tensor. On the other hand, eq. (7-1) below is itself an approximation, 
because the exact field equations should include the Riemann and Ricci terms created by the 
EM field and should be coupled to the gravitational de Rham Laplacian. The formulation 
presented here is based on the hypothesis of test fields (no backreaction) and the goal of 
this analysis is to recover the existing perturbation theory, including the work of Teukolsky, 
Fackerell and Ipser and many other authors. 

The sourcefree Maxwell equations 8F = = dF for the Maxwell 2-form F (representable 
as 4 complex equations for the associated NP scalars (f) , 0i, 2 , namely eqs. (330)-(333) of 
Chapter 1 of 37 ^) immediately imply the vanishing of the de Rham Laplacian of the Maxwell 
2-form. In a vacuum spacetime, this takes the form 

(AdR)^U = F ab \ c \ c + C abcd F cd = , (7-1) 
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Expressing these equations in the Newman-Penrose formalism, one has the following wave 
equations. 



[Z\(dR)-F]i3 = (wave equation for <fio) 

Studying the component "13" of (7-1) we get the equation for fa: 

[DA + AD- 5*5 - 55* + (-5 7 - 7* + p + p*)D + (-3e + e* - p* - p)A 
+ (-(3* - 7t + 5a + t*)5 + (-7T* + r + 3(3 + a*)5* + 2(5* (3 + 5a - L> 7 - Ae 
+ 7*e - 7e* - T7r + 7 p* + ir*a + 4 7 e - (3t* + pp + (3(3* - ep* - aa* - Aa(3 
+ p 7 - kv + (3-k + crA - to - pe - ^ 2 ]0o + [4(r£> + kA - po~ - <x<T) 
+ 2(£>r + - 5p - 5*a + e*r - pvr* - er + err* - «/y* - rp* - cr/?* + pa* 
+ pK + p*K + (3p-a-K- 3«7 + 3ao- + 2^)]^ + [Akt - Aap - 2%}(f) 2 = . 

First we use the commutation rules (2-7); the next step is to divide the resulting equation 
by 2 and to use the two Maxwell equations: 

(D - 2p)<j> 1 - (5* + 7T - 2a)0 o + /t0 2 = 0, (7-3) 
(5 - 2r)0! -(A + p- 2 7 )0 O + ^2 = 0, (7-4) 

to eliminate two of the directional derivatives of <f>\. Then using four Ricci identities: 
[l/2(#i2i2 - ^3412)] : Ae = D 1 - a(r + n*) - (3(r* + vr) + 7 (e + e*) 

+ c(7 + 7*) - ttt + i//c - $2, ^ 
[1/2(^1234-^3434)]: = 5a - (pp - Ac-) - aa* - (3(3* + 2a(3 

- 7(P - P*) - e(p - P*) + #2, 
together with the two in (6-13), rewriting the Ricci identity [-R2431] as in (6-11), multiplying 
it by 0o an d adding the resulting quantity (which is actually zero) to the wave equation one 
obtains: 

[(D — e + e* — 2p — p*)(A + p - 2 7 ) - (5 - (3 - a* - 2r + n*)(5* + vr - 2a)]0 o 

+ (tr A - Kv)<f> + 2[kA - a5* + Ak- 5*a - - a(3* + np* + or* + 3aa (7-6) 

-3«7 + 2^i]0i -!^o02 = 0, 

where in the first line it is easy to recognize the operator found by Teukolsky for the elec- 
tromagnetic perturbations of a type D metric. One can also verify that (7-6) is the electro- 
magnetic equation of Stewart and Walker in the GHP formalism: 

[pp' - d&' - p'p - (2p + p)p' + t'& + (2r + f')& + 2pp' - 2tt' + ^ 2 ]0 O 

+ [2f>'« - 25V - 2re(p' - 2p') + 2<r(f - 2r') + 4^]^ + [2k<3 - 2a\> (7-7) 

+ 2pcr - 2kt - %}<f) 2 = 0. 
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In fact "translating" (7-7) into the standard Newman-Penrose form one gets: 



[(£> - e + e*)(A -2 1 )-{8-(3- a*)(5* - 2a) + /i(D - 2e) 
- (2p + p*)(A - 2 7 ) - n(5 - 2(3) + (2r - tt*)(5* - 2a) - 2pp + 2ttt + ^ 2 ]0o 
+ 2[(Z\ - 37 - 7*)k - (5* -3a + (3*)a + n(p* - 2/i) + <t(t* + 2tt) + 2lPi]^i 
+ [2k(5 + 2(3) - 2o[D + 2e) + 2p<x - 2kt - ^ o ]02 = . 



(7-8) 



Using the following two Maxwell equations: 



(D - p + 2e)0 2 - (5* + 27r)0! + A0 O = 0, 
(5 - r + 2/5)0 2 - (Z\ + 2p)0 x + z/0 o = 0, 



(7-9) 



in the previous equation and adding the Ricci identity [-R2431] (see eq. 6-11) multiplied by 
0o, one finds identically eq. (7-6). 

(wave equation for 2y ) 



l(dRp 
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The equation for 2 can be immediately obtained form eq. (7-6). Using the exchange 
symmetry 1 <-> n, m <-> m* of the Newman-Penrose quantities and the relation (2-18), one 
gets 



- {[(A + 7 - 7* + 2/i + fi*)(D -p + 2e)-(5* + a + [3* + 2ir-T*)(5-T + 2/3)]0 2 
+ (aX - kv)4>2 + 2[X5 - vD + 5X - Du - ue* - Xa* + up* + Atf* + 3/3 A (7-10) 
-3z/e + 2^ 3 ]0 1 -^ 4 o } = . 

The overall minus factor in our case can be eliminated but this form is convenient for further 
developments in the presence of sources. The GHP analogue of this equation is 



[p'p - &'d - pp' - (2p' + p')p + rd' + (2r' + f )& + 2p'p - 2r'r + ^ 2 ]02 

+ [2f>«' - 25a' - 2«'(p - 2p) + 2a' (f' - 2r) + 4^ 3 ]0i + [2«/d' - 2</p' (7-11) 

+ 2pV - 2kV - !^ 4 ] 0o = 0. 



[AdR)^1i2 - [A {dR )F]u = 



(wave equation for <pi) 
To obtain the equation for 0i one has to subtract the components "12" and "34" of (7-1) 
in order to eliminate all the 0* terms. The resulting wave equation, divided by a factor of 
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2, is: 

[DA + AD- 5*5 - 55* + (-7 - 7* + p + p*)D + (e + e* - p* - p)A 

+ {-(3* - vr + a + t*)5 + (-7r* + r - /3 + + 4(pp + Act - ttt - i/k + ^ 2 )]0i 

\-2vD - 2txA + 2A5 + 2p5* - Dv - An + 5X + 5* p - Xa* + pv + tt 7 * - pr* 

(7-12) 

- rA - 3ap + Att* - vrp* + p/3* - Xf3 - ue* + up* + 377r + i/e - 2^ 3 ]0 O 

+ [2tD + 2/tZ\ - 2p5 - 2a5* + Dr + Ak - 5p- 5*a - an - rp* + 7/c - K7* 

+ 3re + rV - 3p/3 + p*/t - atr + re* - a/?* + pa* + pK - pvr* - 2$i]0 2 = 0. 

The procedure to simplify the previous relations is the usual one. The starting point is 
to rewrite the second order part by using (6-6); dividing the resulting equation by 2 and 
using all four NP Maxwell equations, some of the directional derivatives of fa and fa are 
eliminated in favor of the derivatives of fa. Now, using the Ricci identies: (6-10), (6-13) and 

[^2443] : SX - 5*p = v{p - p*) + vr(p - p*) + p(a + /?*) + X(a* - 3/?) - V 3 
[-R2421] : J Dz/-Z\7r = / u(7r + r*) + A(7r* + r) + 7r(7-7*)-z/(3e + e*) + ^ 3 ; 

adding the Ricci identity (6-11) multiplied by 2fa, we finally obtain 

[(D + e + e* - p - p*)(A + 2p) - (5 + (3 - a* - r + n*)(5* + 2ir)]fa 

+ \-vD + X5 - Dv + 5X- Xa* - eV + pi/ - r A + A/3 - z/e + z/p* + n*X]fa (7-14) 

+ — a5* + Z\/t — 5*(X — (7/3* — K7* + p/t — <77T + /tp* + (7T* + «(J — K^]fa = . 

It is easily seen that in the (exact) eq. (7-14) the fa part is exactly the Newman- Penrose 
translation of the GHP equation which describes the fa electromagnetic perturbation in a 
vacuum type D metric present in the paper of Lun and Fernandez 42 ) : 

[(P - P ~ P)(¥ + 2p)-(d-r + n)(d' + 2n)]fa = , (7-15) 

which, after a tetrad transformation which sets e = 0, becomes the well known Fackerell and 
Ipser equation 43 ^. 

It is important to get the general vacuum equation for fa in the GHP formalism, which 
can be quite easily obtanined by following the Stewart and Walker prescriptions. We only 
give the final result: 

[(P - P ~ PW - V) - (d - r - ?){& - 2r')]fa 

+ [k' P + K - a'O - da' + a'r + a'f' - K 'p - n'plfa (7-16) 
+ \k\>' + p'/t — ad' — d'a + fa + r'a — p'n — p 1 n\fa = . 
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Translating (7-16) into standard Newman- Penrose formalism one gets (7-14). 

At this point the (first order) perturbation analysis can be performed. In the type D 
case, recalling that for vacuum spacetimes one has to use the test field approximation (no 
changes in the vacuum background means no changes in the tetrad vectors, spin coefficients 
and Weyl scalars), it is straightforward to show the coincidence of the linearized form of 
(7-6), (7-10) and (7-14) with the Teukolsky, Fackerell-Ipser and Lun-Fernandez equations. In 
a more general vacuum geometry the decoupling is not present and one obviously recovers 
the Stewart and Walker theorems. 

§8. Concluding remarks 

In this article we have presented a new form of the Teukolsky Master Equation and a 
new point of view for perturbations in General Relativity, starting from exact wave equa- 
tions in curved vacuum spacetimes valid for the Riemann and Maxwell tensors themselves 
and involving both the generalized and the ordinary de Rham Laplacian respectively. Once 
perturbed these equations immediately give the standard perturbation theory of Teukolsky, 
Stewart and Walker, Lousto and Campanelli, explaining also the true origin and the geo- 
metrical meaning of the same Teukolsky Master Equation. We have discussed our exact 
wave equations both in the Newman-Penrose and Geroch-Held-Penrose formalisms, collect- 
ing most of results scattered in literature in a series of papers and books and over a period 
of more than thirty years. We have also presented a systematic study of perturbation theory 
at every order, which recovers and extends existing results on this subject. Furthermore, our 
approach can be directly generalized to nonvacuum spacetimes: 1) The de Rham Laplacian 
of the Riemann tensor then involves (in the exact regime) the Ricci tensor and curvature 
scalar terms, in turn to be re-expressed in terms of energy-momentum tensor of the sources 
or directly coupled to the de Rham equations for the other corresponding physical fields: 
Maxwell, Yang Mills, Dirac, Rarita-Schwinger 44 ) ; 2) the perturbation regime can be devel- 
oped straightforwardly, recovering, for instance, once our general results are specialized to 
a vacuum type D background, the "matter perturbed sources" appearing in the works of 
Teukolsky and of Lousto and Campanelli. 

Regarding the tetrad and gauge invariance of the theory, the procedure that one must 
follow is the standard one, well discussed in the literature by Bruni, Matarrese, Mollerach 
and Sonego 45 )' 46 ) and by 35 ) and we will adapt it to this formalism in a forthcoming paper. 

We conclude pointing out that our approach could have useful consequences for the 
still open problem of the electromagnetic and gravitational coupled perturbations of the 
Kerr- Newman spacetime 47 ^ " 50 ) as well as for the more general problem of propagation of 
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gravitational waves on any given background. 
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